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Abstract: Higgs boson and massive-graviton productions in association with two jets via 
vector-boson fusion (VBF) processes and their decays into a vector-boson pair at hadron 
colUders are studied. They include scalar and tensor boson production processes via weak- 
boson fusion in quark-quark collisions, gluon fusion in quark-quark (qq), quark-gluon (qg) 
and gluon-gluon (gg) collisions, as well as their decays into a pair of weak bosons or 
virtual gluons which subsequently decay into ii, qq or gg. We give the helicity amplitudes 
explicitly for all the VBF subprocesses, and show that the VBF amplitudes dominate the 
exact matrix elements not only for the weak-boson fusion processes but also for all the 
gluon fusion processes when appropriate selection cuts are applied, such as a large rapidity 
separation between two jets and a slicing cut for the transverse momenta of the jets. We 
also show that our off-shell vector-boson current amplitudes reduce to the standard quark 
and gluon splitting amplitudes with appropriate gluon-polarization phases in the collinear 
limit. Nontrivial azimuthal angle correlations of the jets in the production and in the decay 
of massive spin-0 and -2 bosons are manifestly expressed as the quantum interference among 
different helicity states of the intermediate vector-bosons. Those correlations reflect the 
spin and the CP nature of the Higgs bosons and the massive gravitons. 
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1. Introduction 

Angular correlation of the two accompanying jets in Higgs boson productions at the CERN 
Large Hadron Collider (LHC) has been known as a potential tool to study its spin and CP 
nature, in the weak-boson fusion (WBF) qq qqH processes [1], and in the gluon fusion 
(GF) plus dijet production processes [2-4], qq — > qqH, qg qgH and gg — > ggH. In these 
reactions, the tensor structure of the Higgs coupling to weak bosons or gluons gives rise 
to the azimuthal angle correlation of the tagging jets; the WBF processes give flat ^(pjj 
distribution, while the GF processes produce a distinct dip around Acpjj = tt/2. On the 
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other hand, in the case of a CP-odd Higgs boson, the azimuthal distribution is strongly 
enhanced around A0jj = 7r/2 in both the WBF and GF processes [1, 3, 4].^ 

So far, many studies on the azimuthal correlations in the Higgs + 2-jet events have been 
performed with higher-order QCD and electroweak corrections [8-13] including parton- 
shower effects [14, 15]. The present consensus seems to be that the azimuthal angle corre- 
lations predicted in the leading order may survive even after higher-order corrections are 
applied [13, 15]. It must be pointed out here that, due to the complicated matrix elements, 
it is not completely clear why the tensor structure of the couplings leads to such distinct 
azimuthal distributions, even though naive explanations have been presented [1, 2]. Since, 
in general, azimuthal angle dependence should be understood in terms of the quantum 
interference phases of the amplitudes with spin-full particles propagating along the polar 
axis, it may be valuable to reformulate the amplitudes for the Higgs production with two 
jets in such a way that their phases are shown explicitly. 

As another interest to study the azimuthal angle correlation of the jets at the LHC, 
we attempt to apply it to other heavy particle productions. Here, we especially focus 
on massive-graviton productions in the localized gravity model of Randall and Sundrum 
(RS) [16], which has drawn a lot of attention in recent years because it brings a new 
solution to the hierarchy problem through an exponentially suppressed warp factor in a 5- 
dimensional non-factorizablc geometry. Several phcnomcnological studies have been made 
on the Drell-Yan process for RS graviton resonances for its discovery and the determination 
of its spin-2 nature [17-23], including direct searches at the Tevatron [24], as well as the 
graviton -|- 1-jet productions [25]. Meanwhile massive-graviton productions in association 
with two jets may also have a great potential to scrutinize its properties as in the Higgs 
boson case. We note that the graviton + 2-jet productions in the large extra dimensions 
model [26] have recently been studied in ref. [27]. 

In this article, more generally, we study productions of a heavy color-singlet particle 
{X) in association with two jets via vector-boson fusion (VBF) processes at hadron collid- 
ers, pp jjX, which include WBF processes in quark-quark collisions and GF processes 
in quark-quark, quark-gluon and gluon-gluon collisions.^ In particular, the reactions, 

qq^qqX, qg ^ qgX, gg ^ ggX, (1.1) 

are studied comprehensively as the leading-order subprocesses that lead to X + 2-jet 
events via VBF. In order to discern the phases of the amplitudes, we present the hclicity 
amplitudes explicitly for all the VBF subprocesses at the tree level in terms of the specific 
kinematical variables, where the colliding vector-bosons have momenta back-to-back along 
the polar axis. Although the VBF amplitudes are valid only when the virtuality of the 
intermediate vector-bosons is smaller than their energies, as we will see later, they can 
dominate the exact matrix elements when appropriate selection cuts to the final states are 

Azimuthal correlations in diffractive processes have also been discussed in refs. [5-7]. 

^Weak-boson fusion is sometimes referred to as VBF. In this paper, however, we refer to the fusion 
processes of all the standard model vector-bosons {W, Z, 7, g) as VBF, including WBF and GF processes. 
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applied, such as a large rapidity separation between two tagging jets and a slicing cut for 
the transverse momenta of the jets. We also show that our off-shell vector-boson current 
amplitudes reduce to the standard quark and gluon splitting amplitudes with appropriate 
phases in the collinear limit. 

As for the produced heavy particles, we study neutral CP-even and CP-odd Higgs 
bosons and RS massive gravitons, and show that nontrivial azimuthal angle correlations 
of the jets in the production of massive spin-0 and -2 bosons are manifestly expressed 
as the quantum interference among different helicity states of the intermediate vector- 
bosons. Those correlation reflects the spin and the CP nature of the Higgs bosons and 
the massive gravitons. We do not consider massive spin-1 particles because the Landau- 
Yang theorem forbids production of a color-singlet spin-1 particle in fusions of two on-shell 
photons or gluons [28], and our approximation fails when their virtuality is large enough 
to give significant amplitudes. 

Besides jet angular correlations in the production processes, the decay distributions 
and correlations of heavy particles are also promising tools to determine their properties, 
and extensive studies have been made especially for the Higgs bosons, e.g. — > ZZ — > 
{££){££) [29-33]; see also review papers [34,35] and references therein. The above decay 
process is related by crossing symmetry to the WBF Higgs production process, qq qqH, 
and hence it may be useful to compare the production correlations with the decay corre- 
lations. Therefore, we also consider X decays into a pair of weak-bosons or gluons which 
subsequently decay into U, qq or gg, and present the helicity amplitudes and the azimuthal 
angle correlations of the jets (and/or leptons), by comparing with those in the production 
processes. 

The article is organized as follows. In section |2|, we introduce the formalism of the 
helicity amplitudes and the density matrices for the X production with two jets via VBF 
and the X decay into a vector-boson pair. In section ^ we define kinematical variables 
relevant to our analysis for the production and decay processes. In section ^, we present 
all the helicity amplitudes explicitly for the off-shell vector-boson currents and the off- 
shell VBF processes. We also discuss the relation between our current amplitudes and 
the standard parton splitting amplitudes. In section |5|, we demonstrate that the VBF 
amplitudes dominate the exact matrix elements when appropriate selection cuts to the 
final states are applied, and then discuss azimuthal angle correlations of the jets in the 
Higgs boson and massive-graviton productions. We also consider the decay correlations of 
the heavy particles. Finally section ^ summarizes our findings. 

We include three appendices. Appendix^ gives the wavefunction and the vertices for 
a spin-2 particle. In appendix^ we show the relation between wavefunctions and Wigner's 
d functions. Appendix ^ presents the angular distributions for the massive-graviton decays, 
G^VV ^{ff){ff). 

2. Helicity formalism 

In this section, we give the helicity amplitude formulae and the density matrix formalism 
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(a) (b) (c) 



Figure 1: Feynman diagrams for the VBF subprocesses: (a) qq — > qqX, (b) qg qgX, and 
(c) 99 99^- 



for heavy particle {X) productions in association with two jets via VBF processes, and also 
those for its decay into a (virtual) vector-boson pair which subsequently decay into ii, qq 
or gg. 

X + 2-jet productions via VBF at hadron colliders, pp jjX, can proceed through 
the subprocesses: 

(2.1a) 
(2.1b) 
(2.1c) 

where V* is a t-channel intermediate vector-boson and q stands for a quark or antiquark of 
any flavors. A representative Feynman diagram for each subprocess is shown in fig. |l|(a), 
(b), and (c), respectively, for the subprocess ( 2.1a| ), (|2.1b| ), and ( 2.1c ). Each subprocess 



qq - 


qqV*V* - 


-> qqX 


{V 


= W, Z, 7, g), 


Q9 - 


qgV*V* - 


-qgX 


{V 


= 9), 


99 - 


ggV*V* - 


-ggx 


{V 


= 9), 



receives contributions not only from the above VBF diagram but also from all the other 
diagrams of the same order, in order to make the gauge- invariant physical amplitudes. In 
this section, however, we consider only the VBF diagram. As we shall see later in section |5|, 
after applying appropriate kinematical selection cuts, the VBF contribution can dominate 
the exact matrix elements. 

Let us first define a common set of kinematical variables for the VBF subprocesses 
( |2.1| ) generically as 

ai{ki,ai) + 02(^2, 0-2) a-siks^a-i) + 04(^4, (74) + V{{qi, Ai) + y2*(g'2, A2) 

^ 03(^3, t^3) + a4{k4, ^4) + X{P, A), (2.2) 

where ai^... ^4 stand for quarks or gluons (or even leptons in case of lepton-lepton or lepton- 
hadron collisions), and the four- momentum and the helicity of each particle are shown 
in parentheses; see also fig. |2|(a). The parton helicities take the values o"j/2 for quarks or 
antiquarks and ai for gluons with = ±1, while the helicities of the off-shell vector-bosons 
take Xi = ±1,0. 
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Figure 2: Schematic view of the subprocesses for (a) the X production with 2 jets via VBF, and 
(b) the X decay to 4 jets via a vector-boson pair. The four-momentum and the hchcity of each 
particle are shown. The solid lines show either fermions or gluons. 



The helicity amplitudes for the VBF processes (^) can generally be expressed as 

"^'^ X ^;:,,('?i)<?,,fe)r^^^,V,(gi,g2;A)*, (2.3) 

where Ji^r „ „ and J(J „ „ are the external fermion or gluon currents, and the vector-boson 
propagators are 

' ' - v. + ^) Dv. (qf) for V, = W, Z, 

^V. ' (2.4) 



-9ti'fi Dv, {qf) for Vi = -f, g 



with the propagator factor Dv{q^) = {q^ — my + imyVy)^^ . Note that we choose the 
unitary-gauge propagator for the massive vector-bosons and the Feynman-gauge one for 
the massless vector-bosons (my. = 0). The XV V vertex is expressed generically as ^^vly^-, 



whose explicit forms are given in section 4.3 



Using the completeness relation for space-like vector-bosons (g? < 0) 

-5m'/. + ^^= E (-l)''+'v('/-A,)*e^(ft,A,), (2.5) 



A,=±,0 



and neglecting the terms which vanish due to current conservation 

9iM'^y,aiai+2(^»'^i+2;^i'^i+2) = 0, (2.6) 



the VBF helicity amplitudes ( p.3| ) can be rewritten as the product of the two incoming 
current (/ fV* or g ^ gV*) amplitudes and the off-shell VBF X production {y*V* 
X) amplitudes summed over the polarization of the intermediate vector-bosons 

V\ 2 Al 2 



- 5 - 



where 



Aside from the summation of Vi_2, the VBF amphtudes ( |2.7| ) are generally the coherent 
sum of the nine amplitudes which have the different helicity combinations of the colliding 



vector-bosons. The explicit forms of the amplitudes (2^) and ( p.9| ) will be given in sec 



tion O and |4.3| , respectively. It is worth noting here that the current conservation (^^) , 
which ensures that the propagating vector-bosons have only the three vector-boson com- 
ponents, plays an essential role in deriving the above expressions, and that it is valid not 
only for currents made of massless fermions but also for those made of on-shell gluons. 

In fig. H, we notice that the diagrams for (a) the X production with 2 jets via VBF 
and those for (b) the X decay into 4 jets via a vector-boson pair have identical topology, 
even though the intermediate vector-bosons are space-like for the production while they 
are time-like for the decay. They are related with each other by the crossing symmetry. 
Therefore, it may be useful to study the decay angular distributions and correlations in all 
the channels simultaneously. 

We consider X decays into a (virtual) vector-boson pair which subsequently decay into 



U, qq, or gg, similar to the Xjj productions via the VBF processes (2.1), 

X ^V'^*'>V'^*'> ^{£e/qq){ee/qq) {V = W,Z,^,g), (2.10a) 
X^V*V* ^{qq){gg) {V = g), (2.10b) 

X^V*V* -^{gg){gg) {V = g), (2.10c) 

and define a common set of kinematical variables as 

x(p',A')-yr(9'i,A;)+yr('?2,A'2) 

^ a'^{k[,a[) + 4(^:^,4) + a'^{k'^,a'^) + a'^{k'^,a'^), (2.11) 

where the same notations for their momenta and helicities are used as in the production 
processes ( p.2p except for primes ('); see also fig. |2|. 

The decay helicity amplitudes can be expressed in the same way as the production 
amplitudes, while the completeness relation for time-like vector bosons has to be 

-5mV + ^^= E V(9.>A^*6^(gU'J- (2.12) 

We note that no extra sign factor for A'j = is needed for q'^ > 0. The outgoing fermion 
or gluon current iV* — > // or gg) amplitudes and the X VV decay amplitudes, corre- 
sponding to eqs. (|2.8D and (|2.9|) for the production, are given by 



(^Ji' ,)a'(7' - ^M(9i>A-) j'^, , , (/c-,/c-_^2;<7i'0-^^2)' (2-13) 

i^vivi)x'x' = ^'xi^kM^l2, a;)* e,M, A'2)*. (2.14) 

^■^12 12 
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The hehcity amphtudes for the X decay processes ( 2.10[) can now be expressed as 



^%^'..W. = EDvM")DvM")Ei^n^^^^^^^ (2-15) 



^1,2 ''^'1,2 



Similar to the VBF production amplitudes ( |2.7| ), the decay amplitudes ( p. 15 ) are generally 



the sum of the nine amplitudes. We note our convention that the sign of the current am- 
plitudes with the longitudinal polarization is different between the incoming and outgoing 
current amplitudes, (^) and ( 2.13| ), reflecting the difference between ( |2.5| ) and (2.12) for 
the spin-1 completeness relation. 

For the sake of our later discussions, we give the complete amplitudes and the squared 
matrix elements for the subprocesses of the X plus n-jet production and its subsequent 
decay into n' jets. In the vicinity of the X resonance pole, the full amplitudes can be 
factorized into the X production amplitudes and its decay amplitudes, summed over the 
helicity of the X resonance. In the X rest frame, if we use the X polarization along the 
z-axis {Jz = A) to express the production amplitudes and that of the decaying direction 
{Jz' = A') to express the decay amplitudes, the full amplitudes can be expressed as 

^^....,..2;<...., = EK...,..2^^(^')<v(e)^l;^; (2.16) 

A,A' ' ' 

where Dx{P^) = {P"^ — + iMT)~^ times the Wigner's d function dj^y^,{Q) gives the 
X propagator in the X rest frame when the X spin is J and the initial and the final 
quantization axes have the opening angle 0. For instance, for the n = case the initial X 
polarization is A = ui — <T2, while the final polarization is A' = ct'i — ct^ for n' = 2. We derive 
the above expression explicitly for the J = 2 case in appendix ^ The amplitudes for the 
scalar particle production and its decay processes have no dependence, i.e. (i^^?(0) = 

dOo(e) = 1. 

It is straightforward to obtain the polarization-summed squared matrix elements of 
the full production plus decay amplitudes (|2.16| ), 

|2 



= \Dx{P^'Y.Y.'^X-AA^)di:y{Q)Vfr,, (2.17) 



A,A' A,A' 



in terms of the production density matrix and the decay density matrix T^yy', 



^A^A'= E ^% X- (2-18b) 

' * I,'-- ,n' I,--- ,n' 



Although eq. ( |2.17 ) applies only for parton-level subprocesses, one can easily generalize it to 



mixed case and apply it for any processes, including summation over different subprocesses 
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and a product of the relevant parton distribution functions. We note that the density 
matrices together with the d functions have all information on the angular distributions of 
the final states. 

For the VBF processes (2.1) and for the X decays into a pair of vector bosons ( 2.1C| ), 



the density matrices can be expressed as 

'Pxx = E \DvMl)Dv,{qlfY.Y.^^^^ - ^2 - X)6{\i - A2 - X)V^,f^, (2.19a) 



^1,2 



^1,2 -^1,2 -^'1,2 



in terms of the colliding or the decaying vector-boson polarization amplitudes, (|2.7| ) or 

^it= E [^l'i.3^2^2.4-^xj;a2][^1^1.3^2^2.4-^xiAj*' (2-20a) 

3. Kinematics 

In this section, we define kinematical variables for the production of the heavy particle 



X and its decay (2.11). The angular configuration of the particles defined below is 
summarized in fig. ^. It must be noted again that the same notations except for primes (') 
are used both for the production and the decay. 

First of all, we note the relations between the momenta by the momentum conservation. 

qi + q2 = P = P' = q'i + q2 (3.1) 

with Qi = ki — ki^2 and q'i = k'^ + k'-^2- On-shell conditions for the external particles and 
the partonic center-of-mass (CM) energy are 

kf = kf^^ = kf = fc-+2 = 0, p2 = P'^ = M^, and s = {ki + ki)^, (3.2) 

while the intermediate vector-bosons are off-shell; 

qf < 0, qf > 0. (3.3) 



We would like to express the VV — > X production amplitudes (|2.9D and the X VV decay 



amplitudes (2.14) in the X rest frame where the colliding or decaying (virtual) vector-boson 
momenta are chosen as the polar axis, along which their helicities are defined. Because the 
helicities are invariant under the boost along the polar axis, we evaluate the vector-boson 
emission amplitudes in the Breit frames, qi^ = \J —q\ or q^^ = — \/— and the vector- 
boson decay amplitudes in their rest frames, (f^ — \J q'l^ q'2^ - 

For the X + 2-jet production process via VBF (^.2|), we parametrize the momenta of 
the quarks, antiquarks or gluons of the incoming currents in the Breit frame as follows: 
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z 



Figure 3: The momentum and angular configuration of the particles in the q\ and Breit frame, 
(I) and (II), and the VBF frame (III) for the production; in the and cl^ rest frame, (I') and (II'), 
and the X rest frame (III') for the decay. 

(I) the q\ Breit frame 

^1 = ^1 - ^3 = (0' 0' 0' Qi)' 

^1 ~ 2ros6'i sin 6*1 COS , sin ^1 sin , cos6'i), 

^3 ~ 2cos6>i (-*^' sin 6*1 COS (/>i , sin 6*1 sin , — cos^i), (3.4) 

where Q\ = \J -q{, < 6*1 < 7r/2 and < 0i < 27r. 
(II) the (72 Breit frame 

= kl^-k^ = (0, 0, 0, -Q2), 
^2 = ~ 2ms 92 sin ^2 COS 02, sin 6*2 sin (/)2 , cos 6*2), 
^4 = ~ 2cos6>2 sin 6*2 cos 02, sin 6*2 sin 02, - cos 6*2), (3.5) 

where (^2 = V ^/^ < O2 < tt and < 02 < 27r. 

The momenta of the t-channel intermediate vector-bosons are chosen along the z-axis, 
and the coordinate system where the vector boson has a positive (negative) momentum 
along the z-axis is labeled as 1 (2). In practice, we always denote by ki the parton 
momentum in the colliding proton with pz = \/i/2, while k2 denotes the parton momentum 
in the other proton (or anti-proton) with pz = —y/s/2, in the laboratory frame. This agrees 
with the above definition whenever the VBF approximation to the full amplitude is valid. 
The polar angles 0i and 62 are measured from the common positive z-axis; see also fig. |3|. 
By a boost along the z-axis, each Breit frame can be transformed to the rest frame of the 
heavy particle, referred to as the VBF frame: 
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(Ill) the VBF frame 

q^i+q^2=P^ = P'" = q'i" + = (M, 0, 0, 0), 
= M(i_Ql_Q|,o, 0, 
g^ = f (1-^,0, 0, -p), 
g;^ = f + /3'smG, 0, /3'cose), 

^/ = M (1 ^ -/?' sin e, 0, -P' cos e) , (3.6) 

where Q[ = /? = ^( - -g, and /?' = p{%, %) with /3(a, 6) = (1 + + 6^ - 

2a — 2b — 2a6)^/^, and is the angle between the production axis (z-axis) and the decay 
axis (z'-axis). The boost factor along the z-axis from each Breit frame to the VBF frame 
is, respectively, 

/3i = (l-%^)//5, (3.7a) 
/32 = -(l-%^)//3. (3.7b) 

We note that, when the produced particle X decays into visible particles and its mo- 
mentum is reconstructed together with those of the two tagging jets, the Breit frames can 
in principle be reconstructed in experiments. 

For the X decays into 4 jets via a (virtual) vector-boson pair ( ^.ll| ), the momenta of 
the time-like vector-bosons are measured along the z'-axis in the X rest frame: 

(III') the X rest frame 

< + q^2 =P^ = P'^ = + 92^ = (M, 0, 0, 0), 
< = f (1 - -/3sine, 0, /3cose), 

12 = f (1 - /9sine, 0, -/3cose), 

gr = f (l + ^^g^,0, 0,/3'), 

9^'' = f(l + ^%^, 0,0, -/?'), (3.8) 

which is obtained from the VBF frame ( ^.(j ) by the rotation with the angle (—0) about 
the y-axSs. 

For the outgoing fermions and gluons, their momenta are parametrized in the rest 
frames of the time-like vector-bosons as follows; 

(I') the q'^ rest frame: 

g'/ = A.'/ + A:^'^ = (Q;,0, 0, 0), 

k'l^ = %-(!, sin 0'^ cos 0']^ , sin 0^ sin ^'^^ , cos0^), 

k'/ = ^(1, -sin6''iCOS(/);, -sin6'isin^;, -cos6''i), (3.9) 
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(II') the I'sst frame: 

q'/ = k'/ + k'/ = iQ'^, 0, 0, 0), 

k'2'^ = ^(1, sin ^2 cos 02, sin02sin02i cos02)) 

y^fi _ ^(^i^ — sin 6*2 cos 02, — sin 6*2 sin 02 , — cos6'2), (3.10) 

where < 9'^ < ir and < 0^ < 27r. In the X rest frame, the q'l momentum is chosen along 
the z'-axis, and the (72 momentmn along the negative z'-axis. The polar (2;'-)axis and the 
y-axis normal to the scattering plane are chosen common to all the three frames for the 
decay chain, III', I' and II', which are related with each other by a boost along the z'-axis; 
see also fig. ^. For instance, both cos6[ = 1 in ( ^!g| ) and cos02 = 1 in ( 3.10 ) denote the 



momentum along the z'-axis direction, and hence the a'l momentum (k'l) is along the V( 
momentum while the 02 momentum (^2) is anti-parallel to the V2 momentum in the X 
rest frame. The boost factor along the z'-axis from each vector-boson rest frame to the X 
rest frame is, respectively. 



/?;=/?'/(! + (3.11a) 
= -/?'/ (l + ^%#i^). (3.11b) 



4. Helicity amplitudes 

As we have shown in section |2|, the VBF helicity amplitudes ( p.7| ) can be expressed by 
the product of the two incoming current amplitudes and the VV X fusion amplitudes. 
Similarly, the decay helicity amplitudes ( 2.15| ) can be given by the product of the X VV 



decay amplitudes and the two outgoing current amplitudes. In this section, using the helic- 
ity amplitude technique [36] and the kinematical variables defined in the previous section, 
we present all the helicity amplitudes explicitly for the fermion currents, the gluon cur- 
rents, and the off-shell VBF vertices, respectively. We also discuss the relation between the 
off-shell vector-boson current amplitudes and the standard parton splitting amplitudes [37] . 

4.1 Current amplitudes 

Let us start with the helicity amplitudes for the incoming fermion currents (/ — > fV*) 
in the X production process, {J^^ai+2)a'a+2 ™ "^^^ fermion and antifermion 

currents are given by 

Jyjf,{ki,ki+2;o-i,ai+2) = g^/^' Uf>{ki+2,<yi+2)l^ Uf{ki,ai), (4.1) 

Jyfp{ki,ki+2]cri,ai+2) = g^l! Vf{ki,ai)j''vf{ki+2,cri+2)- (4.2) 
Non-vanishing couplings in the standard model (SM) are 
9± — eQf, g± — Qst , 

gZ^ff ^ -gzQf sin^ Ow, g^^^ = gzlTsf - Qf sin^ 9w], 

Wu.,, ^ ^Wd,u.y ^ (4.3) 
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[cos 6*1 — 2:1/(2 - zi)] 


J+=-{j^- )* 
"-^ '-++ V*^ — / 
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2 cos 61 zi 

1 . a V2(l-^l) 
sill C/l — 

\/2cos6'i 

^ (1 cos 6^1) e*-?^! I~^iei0i 
2 cos 9i zi 






[cose'j 2z[ - 1] 


^iT- = -(^i-+) 

— -'l — V 

— ^1 


- ( 1 + cos 6[ ) 6**^1 4 e*"^! 

-^sin^i ^24(1- 4) 
^(1 - cos e[) e-< (1 - 4) e-^<^'i 




Table 1: The reduced helicity amplitudes for the ofF-sheU vector-boson currents: the incoming 
fermion currents J\^\ajK_f — > /^*) in the Breit frame (top), and the outgoing fermion currents 
J\i) a' (^* ~^ ■^^) vector-boson rest frame (bottom). In the third column the splitting 

amplitudes are also shown in the coUinear limit, where is the energy fraction of the initial 
particle. 



where e = -v/4vra is the magnitude of the electron charge, Qs = y^lvraj is the QCD cou- 
pling constant, e = gwsin9]y = (7^ sin cos j is the SU{3) color matrix, and Vij 
denotes the Cabibbo-Kobayashi-Maskawa (CKM) matrix element. Using the kinematical 
variables defined in the previous section and contracting the above current with the final- 
state polarization vector e^{qi, Xi)* , we can obtain the helicity amplitudes explicitly for all 
the helicity combinations; see eq. ( |2.8D . For our analytical calculations, we use the HELAS 
convention [38] for the spinors. For the virtual vector bosons with space-like momentum, 
qf = —Qf < 0, we define the longitudinal polarization vectors as 

e''iq^, \ = 0) = m'/q^, q, ) (4.4) 

with Qi = \J\{q^Y ~ l^iPl- By choosing the transverse polarization vectors e^(gi, Aj = ±) 
about the qi axis, the identity ( |2.5D holds. 

Here, for notational convenience, we define the reduced current amplitudes ■Ji^\ai^2 ^ 

(^..J t... = ^^'"'^'Q^ (4-5) 

In table |l](top), the reduced helicity amplitudes for the incoming fermion currents, 

(/^^ — > /o-3^j)) are shown in the qi Breit frame ( |3.4[) , or the frame I in fig. ^. The following 

features of the amplitudes are worth noting: (i) The reduced amplitudes for the antiquark 



- 12 - 



currents are the same as those for the quark currents, (ii) Parity transformation gives the 
relation jfi^i.c-a = (J^i-al -(73)* ■ (iii) The quark masses are neglected (eq. (|3^)), 

and hence the helicity-flip amplitudes are zero due to the chirality conservation. 

(iv) The amplitudes ^2 are related to ji by 

Mla,i02,h) = i-l)^' + 'jl^,i\ie2,C^2). (4.6) 

(v) The 1 / cos 9i dependence comes from the common factor of the four- momentum in 
the qi Breit frame (^^), and this gives rise to the enhancement of the amplitudes when 
cos^i approaches zero. This singularity simply reflects the well-known soft singularity in 
the laboratory frame, and will be discussed further in section |4.2| . (vi) The transverse 
currents JT^ have opposite phases with each other in terms of the azimuthal angle, while 
the longitudinal currents do not have the azimuthal angle dependence. 

Next, for comparison, we consider the helicity amplitudes for the outgoing fermion 
currents (V* — > //) in the X decay process, , in eq. ( ^.131) . The fermion 



currents in which the time-like vector-bosons decay into ££ or qq are given by 



The reduced amplitudes for the outgoing currents, defined as in eq. (4.5), 



can be obtained by contracting the current ( [4.7^ with the initial-state polarization vector 

e^((Z^,A'J; see eq. Q). 

In table [^(bottom), the reduced helicity amplitudes for the outgoing fermion currents, 

Jl^} I {Vxi — > / ifi), are shown in the q\ rest frame (^), or the frame I' in fig. 0. The 

amplitudes have similar features to the incoming current amplitudes in table ||(top). It 

"A' - 

is worth noting that not only the s-channel amplitudes JiJ^i (V* ff) but also the 
t-channel ones J^ialasif ~^ f^*) can be expressed by the same J = 1 d functions as 

JilUU ^ fv*)^dl^„^^^.^{e,), (4.9) 

The different points from the incoming current amplitudes are the following: (i) There is 
no 1/ cos 01 factor, (ii) The phases have opposite signs, reflecting the wavefunctions of the 

'■ A' 

incoming and outgoing vector-bosons, (iii) The same-helicity amplitudes „i are zero, 
which correspond to the chirality-flip amplitudes, (iv) The amplitudes J2 are given by 

^4<(^2,0'2) = -^l';^4(^^,'A'2). (4.11) 

It is worth pointing out here that the list of the current amplitudes in table || is useful 
not only for hadron colliders but also for e"'"e~, and 77 colliders. 
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Turning now to the gluon current amphtudes, the incoming gluon currents in the X 
production process, g — > gV* , where V* is a virtual gluon, are given by 

X [ - g'-P^ki + h+2r - + Qir - a'^'i-qi - hf] , (4.12) 

where f"'^'^ is the structure constant of the SU{3) group. Similar to the fermion currents, 
the reduced amplitudes are defined as 

= ^5.r''^Q^ii^W,. (4.13) 

For the polarization vectors of the external gluons in the amplitude ( [4.13 ), we adopt a 
common light-cone gauge, nf = (1, — ^i/|^i|), which satisfy 

rii ■ e{ki,ai) = nt ■ e{ki+2, 0-4+2) = nf = 0, m ■ ki ^ 0, . ki+2 ^ 0. (4.14) 

It should be noted that these gauge-fixing vectors are boost invariant along the current 
momentum directions, in particular between the Breit frames of I and II for i = 1 and 2, 
respectively, and the collision CM frame III. 

By the crossing symmetry, the outgoing gluon currents in the X decay process, V* 
gg, are obtained by making the replacements in eq. ( 4.12 ): ki — > — A:^, ki+2 — *■ ^+2; 
Qi -Qi, and e„ e* . 

In table ^, we present the reduced helicity amplitudes for the incoming gluon currents, 

Ji^iCTa (^CTi ~^ 5(73 ^Ai)i the qi Breit frame (top), and for the outgoing gluon currents, 
* A' 

J\^) ^1 {Vy 9 a' 3 a' )' '^1 ^^^^ frame (bottom). The derivations are straightforward 

as in the fermion case and the amplitudes have similar properties to those for fermions 
in table |l[ However, the results for the gluon currents are more involved. The incoming 
amplitudes ^2 for the helicity-conserved currents ((T2 = <^a) are given by eq. (|4.6|), while 
the outgoing ones J2 for the opposite-helicity currents {a'2 = —(^4) are obtained by 

:^2't(^2,02) = ^l';5(^2,'A2). (4.15) 

Unlike the fermion current amplitudes, some of the helicity- flip amplitudes (ui = —(73) for 
the incoming currents and the same-helicity amplitudes {a[ = Ug) for the outgoing currents 
are nonzero, and the amplitudes ^Tg'^ are given by 

'^2 <^4 <^2 '^4 

Furthermore, in addition to the singularity of the amplitudes at cos 9i = 0, which appears 
also in the incoming fermion amplitudes in table |^(top), the singularity at sin^|'^ = 0, or 
cos 9^^ = 1 also exists; see more discussions in section 14.21. It must be stressed here that 



the phase dependence of the gluonic currents and that of the fermionic currents are very 
similar, or e^^'^'^''^ when the gluon helicity is conserved (g? < 0) or flipped {q'1 > 0) 

as in the fermionic case. These phases lead to the azimuthal angle correlations of the jets, 
as we will show later. 
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-— i -(1 + cos e,f e-*^i i e-^<^i 
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a/2 cos 6*1 \/2zi 

■ / , (1 - COS e^^^ -^^^^== ^'"^ 
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[cose[ 2z[ - l] 


J\^- — ^1-+ 

>-'l++ — •-'1 — 


. (1 +cosf^')2p^<Ai ^1 e*0i 
2sm0'/ ^z[il-z[) 

1 2zi - 1 

= cos Of ^ — 

a/9 a/9 

^ „(1 cosr^;)2p-i<Ai 





Table 2: The same as table |l], but for the gluon currents, g — > gV* (top) and V* — > (bottom), 
where V* is an ofF-sheh ghion. 

4.2 Relation to the splitting amplitudes 

Before turning to the XVV amplitudes, it may be valuable to discuss the off-shell vector- 
boson current amplitudes from a different point of view, parton branching description [37] , 
where the outgoing particles are emitted collinearly. 



To begin with, we consider the incoming current amplitudes (/ — > fV* or 17 — > gV*). 
Let z be the energy fraction of the initial parton that is carried off by the space-like vector- 
boson. In the VBF frame ( |3.6| ), the energy fraction zi = qi/k^ is written in terms of cos 9i 
defined in the qi Breit frame (3^) and the boost factor /?i in eq. ( |3.7a ) as 



^^^JAcos^ (417) 

' k<i i+Acos^r ^ ' 

Taking the /3i = 1 limit, where the space-like vector-boson becomes on-shell and collinear 
with the final parton, we obtain the simple relation between the Breit-frame angle cos 9i 
and the energy fraction zi as 

cos6'i = — — > . (4.18) 

/3i(2-zi) 2-zi ^ ^ 
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In the third column in tables |l|(top) and ^(top), using the above relation ( 4.18| ) in the 



/?! = 1 hmit, or in the collinear hmit, the helicity amplitudes for the incoming fermion and 
gluon currents are rewritten as splitting amplitudes with appropriate phases. These formu- 
lae may give us clear explanation for the origin of the singularities of the amplitudes which 
we encountered in the previous section. We see from table |^(top) for the incoming fermion 
splitting that the amplitudes are enhanced at zi 0, where the vector boson becomes 
soft. On the other hand, for the g gV* splitting in table |2|(top), the enhancements of 
the amplitudes at zi — > and 1 are associated with the soft emissions of the space-like 
gluon and the outgoing jet gluon, respectively. 

Next, we consider the outgoing current amplitudes {V* — > // or gg). Here, we define 
a fraction z' as the energy transferred from the time-like vector-boson to the outgoing 
fermion or gluon. In the X rest frame (p.8|), 

, k'' 1 



,V = TT^(l + /3lcos0'i), (4.19) 



2f3[ 



where cos 6*'^ is defined in the q'l rest frame ( ^^ ) and P'l is the boost factor in eq. ( |3.11a| ). 
In the P'l = 1 limit, where the time-like vector-boson becomes on mass-shell and the two 
outgoing partons are emitted collinearly, we obtain 

cos9[=2z[- ^ 2z[-l. (4.20) 

Pi 

In the third column in tables |l|(bottom) and |2|(bottom), by making the replacement 



( 4.20 ) in the f3[ = 1 limit, we present the splitting amplitudes for the outgoing fermion 
and gluon currents. There is no singularity for the V* — > // splitting amplitudes, while 
the V* — > gg amplitudes have the singularities at z[ = and 1, similar to the space- like 
gluon splitting amplitudes. This is because the singularities are associated only with soft 
gluon emissions. 



Finally, let us confirm that the splitting amplitudes discussed above reproduce the 
standard (unregularized) quark and gluon splitting functions. From table |l|(top), the sum 
of the squared amplitudes for the transversely-polarized vector-bosons (Ai = ±1), aver- 
aged/summed over initial/final parton helicities, leads to the / — > splitting function 
as 

Z Zl 

o^l,3=± 

where the extra factor zi comes from the initial-state flux factor for the space-like branch- 
ing. Similarly, from table |2|(top) we can obtain the gluon splitting function as 

^17 E [I:^i+.3l' + fc.3n =^ + T^ + ^i(1-^i) = ^Vt/i^t(^i)' (4-22) 

t Zl L Zl 

1^1, 3=± 
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where, in addition to the spin averaged factor, we divide by the statistical factor for the 
two identical gluons in the final state. On the other hand, from the V* — > // splitting 
amplitudes in table ||(bottom), which are time-like branching, the Vr — > / splitting function 
can be reproduced, 

I E + I^M^l'] = ^'i' + (1 - -i') = Pf/Vri^i)- (4-23) 



Likewise, from table |2|(bottom) , we can obtain the gluon splitting function as in eq. ( 4.22| ), 



\ E [I^M.^l' + I^M.il']=^v^T/VT(4)- (4.24) 

In the parton branching description, the space-like and time-like vector-bosons, i.e. 
gluons, are almost on mass-shell, and hence their polarization vectors are taken to be 
purely transverse. However, the longitudinal component of the polarization (Ai = 0) also 
exists for the massive vector-bosons. Therefore, in addition to the above standard parton 
splitting functions for the transversely-polarized vector-bosons, we list the functions for 
the longitudinal polarization: 

^4 E l^i°i-3l' = ^^^— ^ = ^v^//(^i) from table l(top), (4.25) 

0'1,3=± 

E I^i°i-3I' = ^^^7^ = ^v^/V't(^i) from table l(top), (4.26) 

E \'^Zwf=Ha-^i) = Pf/vA^i) from table l(bottom), (4.27) 

I E \'^Z[.'f = liH-'^? = Pvr/vM) from table l(bottom). (4.28) 



In fact, many studies have been performed for Higgs boson productions via WBF in the 
equivalent weak-boson approximation, where the t-channel intermediate weak-bosons are 
viewed as partons in the incoming quarks and the above splitting functions -fV^// ^ well 
as Pvt/J considered [39]; see more details in ref. [34] and references therein. 

4.3 OfT-shell VBF amplitudes 

We will now show the final piece, the VV X production and the X — > VV decay 
amplitudes in eqs. (^^) and ( p. 14 ), respectively. In this paper, we consider the productions 



and the decays of massive spin-0 and spin-2 bosons: neutral CP-even and CP-odd Higgs 
bosons {X = H and A), and graviton resonances {X = G). 

For the fusion vertex of a CP-even Higgs boson, we consider both the WBF process 
and the GF process through a top-quark loop in the SM, and their tensor structures 
^xViV2^^^^^'^^^ normalized by the coupling form factors 5'xviV2 ('^i' ^2), are given in table |3|. 
The constant coupling, g^^yy = 2my/v {V = W,Z) with v = 246 GeV, gives the WBF 
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(A) 






H 


(0) 


w,z 




H 


(0) 






A 


(0) 






G 


(±2, ±1,0) 


W,Z,j,g 





Table 3: The XVV vertex T'^y^y^^qi, q2; X) in eq. (p^), normalized by the scalar form factor 
f Js:ViV2 (^i' 92), are defined for CP-even and -odd Higgs bosons {H and A) and massive gravitons 
(G), respectively. For the polarization tensor e°'^{P,X) and the GVV vertex f ^y^^'^^ (gi, (72), see 
appendix |^. 

vertex, while the explicit expression of the form factor gugg{qi,q2) by a triangle-loop is 
given in ref. [2].^ Note that in refs. [1, 3, 8, 10] the same loop-induced vertex structure has 
been considered to study the anomalous couplings between the Higgs and weak bosons. 
We also consider the GF vertex for a CP-odd Higgs boson, defined in table |3|. For light 
Higgs bosons {Mh,a < '^mt), the above Hgg and Agg vertices can be well described by the 
heavy-top effective Lagrangian [40-43] 

^H,A = -\gHggHF^^,F-'^^ + lg^,,^F;,F'^''^^ (4.29) 

where F"-'^^ is the gluon field-strength tensor and F"-'^^ = ^e'^'^P^Fp^ is its dual. The 
coupling constants are given by g^gg = ctsgjj^^/3-Kmt and g^gg = asgy^^i/2'irmt. The same 
tensor structures can be written for the interactions of the CP-even/odd Higgs boson with 
two photons [44] and with a Z-boson and a photon [45]. 

For graviton resonances, we adopt the simplest RS model [16, 17], where only gravitons 
can propagate into the extra dimension, and consider the first excited mode of the Kaluza- 
Klein (KK) gravitons. The low-energy effective interactions with the SM fields are given 

by 

>CG = -^r^"G^„ (4.30) 

where T^'^ is the energy-momentum tensor of the SM fields (see, e.g., ref. [46] for the ex- 
plicit forms) and Gfj_iy is the spin-2 KK graviton. The RS graviton excitations have the 
universal coupling strength of 1/A to the matter and gauge fields, e.g. gQyy = — 1/A (y = 
W, Z, 7, g), where A is the scale parameter of the theory and can be a few TeV. The explicit 
forms of the polarization tensor for a spin-2 graviton and the three-point GVV vertices are 
given in appendix 



^^Strictly speaking, the Hgg coupling tensor is given by two terms as T^i''^ = F^Tf^^^'^ + F^T^^^'^ [2], 
where TJj^^'^'^ is identical with the form in table ^. However, we neglect the second term since it does not 
contribute in the on-shell gluon limit and is not enhanced in the coUinear limit. 
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A (A1A2) 


CP-even CP-odd 
P'(WBF) (loop-induced) A 


(±±) 
(00) 


-1 -^(M2 + g2 + g2) ^l^^^M-^^Qi^Qif-^QiQi 



Table 4: The reduced helicity amplitudes for the CP-even/odd Higgs boson productions via 
oflF-shell vector-boson fusion, Mx\^x^{V^V^^ — > H/A^), in the VBF frame. M is the Higgs boson 
mass, and Qi and Q2 are magnitudes of the four-momentum squared of the vector bosons. 



In table ^ we present the reduced helicity amplitudes for the CP-even/odd Higgs 
boson productions via off-shell vector-bosons, M.x\^\^{y\^V^_^ — > H/A^), where 

{^v,v,)lx, = ^Mife, Ai) 6^,(92, A2) r^^^/^,V.(9i,'?2; A)* 

= 9xv^V2iQi^Q2)Mxx,x,, (4.31) 

in the VBF frame. Since Higgs bosons are spin-0 particles (A = 0), we have the helicity 
selection rule (Ai = A2), and hence only three amplitudes among the nine amplitudes which 
have the different helicity combinations of the colliding vector-bosons can be nonzero. In 
the table, the amplitudes are expressed in terms of the Higgs boson mass M and magnitudes 
of the four- momentum squared of the vector bosons Qi and Q2- The amplitudes for the 
CP-even and CP-odd Higgs bosons via the collisions of the transversely-polarized vector- 
bosons have the relationships, respectively, 

Mh++-Mh~~=0, (4.32) 
Ma+++Ma~- =0- (4.33) 

Moreover, the pseudoscalar Higgs bosons {A) cannot be produced through the longitudinal 
vector-bosons, TWaoo = 0, due to the CP-odd property. We should notice that, in the case 
of Qi,Q2 *C M, the Higgs bosons via WBF are produced mostly through the longitudinally- 
polarized weak-bosons. On the other hand, the loop-induced CP-even (CP-odd) Higgs 
bosons are produced mainly (only) by the transversely-polarized vector-bosons, and the 
magnitude of their amplitudes are almost equal, ~ I-^aIj apart from the overall 

coupling factors. These characters of the XVV amplitudes, together with the phases of 
the current amplitudes, play an important role to develop the distinctive azimuthal angle 
correlations of the jets. 

Similarly, table || shows the reduced helicity amplitudes for the massive-graviton pro- 
ductions via off-shell vector-bosons, ■M.x\^x^iyx^y\^ ^a)' VBF frame, where 
A = Ai — A2 is the tensor helicity along the colliding vector-boson axis (positive z-axis). 
For the spin-2 particle productions, the amplitudes in all the helicity combinations of the 
vector bosons exist, that is, the nine amplitudes are registered in table ^. However, in the 
case of Qi, Q2 the A = ±2 states of the gravitons are dominantly produced for the 
massless vector-boson collisions, when the vector bosons have the opposite-sign transverse 
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A 




±2 


(±t) 


1 1 


(±Uj 


±1 


(OT) 





(±±) 





(00) 



G 

-{M^ + Ql + Ql + 2ml) 



1 



1 

V2MQ1 
1 



[Q?(Af2 + Ql- Ql) - ml{M^ -Ql + QD] 



[{Qj-Ql? + MHQl + Ql-2ml)] 



V6Q1Q2 



4QlQl + 2ml,{M^ + Qi + Ql) 



-^{{M' + Qj + Qir-AQjQl} 



Table 5: The same as table ^, but for the massive-graviton productions, V^_^V^^ G^, where M 
and my is the graviton and the vector-boson mass, respectively. 



polarizations. For the massive vector-boson case, on the other hand, the amplitudes of 
other three states (A = ±1 and 0) are not negligible. 

The X VV decay helicity amplitudes in eq. ( |2.14D , {-M' yiyi)\i y ' obtained from 
the above VV X fusion amplitudes (tables ^ and ^) by making the replacements — > 
—Q'i^, where A'^^ and A2 are the helicities of the decaying vector-bosons and A' = \'i — A2 
is the X helicity along the momentum direction of the vector bosons (z'-axis in fig. |3|). 
We note that for the non-scalar particle decays the dependence of the angle between 
the initial polarization A along the z-axis and the final polarization A' along the z'-axis is 
dictated by a d function and factorized in our convention; see eq. ( 2.16| ). 

At this point it is worth noting that the amplitudes of the Higgs boson decays into 
transversely-polarized weak-bosons {M'ff±±) are comparable to that into longitudinal ones 
(TW'^oo) near the VV threshold, while the longitudinal amplitudes are dominant for large 
Higgs boson masses. This is because one or both of the weak bosons can become on- 
shell and the typical mass scale of the decaying weak-bosons are Qi ~ my. The helicity 
amplitudes of the Higgs boson decays into the on-shell weak-bosons are given by 



M', 



H±± 



-1, M', 



HOO 



1 = 7'2(l + /3'- 



(4.34) 



where /?' = y ^ ~ 4my/M^ and 7' = 1/ yjl — (3'"^ are the velocity and the Lorentz-boost 
factor of the weak bosons, respectively.^ 

It may be also useful to present the helicity amplitudes for the graviton decays into a 
on-shell vector-boson pair; see Table |6|, which can be reduced from the VV G amplitudes 



in table g with the replacements Qj 



-m 



V 



In the heavy graviton- mass limit, or in the 



/?' = 1 limit, only three (two) amplitudes among the nine survive for the decays into massive 
(massless) vector-bosons. 



*See ref. [32] for the amplitudes of the Higgs boson decays into a pair of virtual and real weak-bosons 
with the most general HVV vertices. 
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A' (a;a'2) 


G 


±2 (±T) 
±1 (±0),(0t) 
(±±)/(00) 


-^(l-/3'2)MV--ig(2-/?'2)M2 



Table 6: The reduced helicity amplitudes for the massive-graviton decays into a on-shell vector- 
boson pair, M'xx'x' {Gy ~^ in the graviton rest frame. 



5. Azimuthal angle correlations 

We are now ready to present the helicity amplitudes explicitly for the X productions with 2 
jets via VBF, eq. (^), and the X decays into 4 jets via a vector-boson pair, eq. ( p.l5 ). In 



this section, we demonstrate that the VBF amplitudes dominate the exact matrix elements 
by appropriate selection cuts to the final state [47] , and then show that nontrivial azimuthal 
angle correlations of the jets in the production are manifestly expressed as the quantum 
interference among different helicity states of the intermediate vector-bosons. We also 
discuss the X decay angular correlations. 

5.1 The VBF amplitudes vs. the full amplitudes 

Although we have considered only the VBF diagrams in fig. ||, there are other crossing- 
related diagrams which have to be taken into account. Representative Feynman diagrams 
for each subprocess, including the VBF diagrams, are shown in fig. ^. 

The Higgs bosons are emitted from each of the circle points in the diagrams, (a) qq — > 
qqH/A: The Higgs boson is radiated off the weak-boson propagator, or the gluon (photon) 
propagator through a top-quark triangle- loop. There is only a t-channel VBF diagram 
in the case of the different-flavor initial state, while an additional s-channel diagram or 
a u-channel diagram exists in the identical-flavor case, (b) qg — > qgH/A: In addition 
to the H/Agg coupling along the gluon lines, the process receives the H/Aggg box-loop 
contribution on the gluon three-point vertex. There are in total 8 diagrams in the large-mj 
limit, ^ in which the top-quark loop is replaced by the effective coupling, (c) gg ggH/A: 
Besides the three corresponding gg gg processes in fig. ^(c), the li-channel diagram 
exists. There are 26 diagrams in the heavy-top limit, including one diagram with the 
effective Hgggg vertex induced via a pentagon loop.^ 

The KK gravitons are emitted from both the circle and the square points in the di- 
agrams in fig. ^, due to their universal couplings to the matter and gauge fields. There- 
fore, the graviton productions have many more diagrams even in quark-quark scatterings 
(fig. 1(a)). 

The point which must be investigated here is whether the VBF amplitudes can dom- 
inate the exact matrix elements, in which all the possible diagrams contribute. Our key 

^Away from the heavy-top hmit, the number of diagrams increases because of the ordering of the gluon 
momenta along the top-quark loop [2]. 

®A CP-odd Higgs boson does not have the Agggg vertex due to the anti-symmetric nature of the 
coupling [43]. 
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Figure 4: Representative Feynman diagrams for the subprocesses, (a) qq qqX, (b) qg qgX, 
(c) 99 99X, which contribute to the X + 2-jet production at hadron coUiders, pp — > jjX. The 
Higgs bosons are emitted from each of the circle points in the diagrams, while the KK gravitons 
are emitted from each of the circle and square points. 

observation is that this happens when we select those events which satisfy the character- 
istic kinematical structure of the VBF processes. Due to the t-channel propagators of the 
vector bosons in the VBF amplitudes in eq. ( p.?] ), the Xjj events via VBF are dominantly 
produced when both Qi { = ) and Q2 { = \/—Q2 ) small. In other words, the 

intermediate vector-bosons in the VBF processes tend to carry only a small fraction of the 
initial parton energies. For small Qi and Q2, the initial partons scatter to far forward and 
far backward, and the heavy particle X is produced centrally. Therefore, the two jets have 
the large rapidity gap, which is often used as the so-called WBF cut to enhance the WBF 
Higgs productions. It should be stressed that this kinematical feature is not particular to 
the WBF production processes, but the QCD productions via the f-channel GF processes 
in qq, qg and gg collisions also have the similar kinematical structure. This suggests that 
some kinematical cuts, such as a large rapidity separation between two jets, may select the 
VBF diagrams dominantly among all the possible diagrams. 

Let us demonstrate numerically that the dijet large rapidity separation is an effective 
kinematical cut to select the VBF amplitudes among the full amplitudes. As the minimal 
selection cuts on the final-state partons, we impose the following kinematical constraints 
for the LHC, required by the detector and jet algorithms: 

PT, > 20 GeV, \7]j\< 5, Rjj = ^Aij]^ + AcP^.^ > 0.6, (5.1) 

where pTj and rjj are the transverse momentum and the pseudorapidity of a final-state 
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C^VEF exact 


Arijj > 3 Arijj > 4 A??,,- > 5 


qq ^ qqH/A/G 
qg ^ qgH/A/G 
gg ^ ggH/A/G 


1.00/1.00/1.58 1.00/1.00/1.43 1.00/1.00/1.25 
1.07/1.05/1.30 1.04/1.03/1.18 1.02/1.02/1.11 
1.07/1.06/1.16 1.04/1.04/1.11 1.02/1.02/1.07 



Table 7: Ratio of the VBF contribution to the cross section with the exact matrix elements, 
cvBF /cexact , for cach subprocess at the LHC, after imposing the inclusive cuts (5.1) and the VBF 
cuts ( |5.2| ) with Arjjjnun = 3, 4 and 5. 



parton, respectively, and Rjj describes the separation of the two partons in the plane 
of the pseudorapidity and the azimuthal angle. Moreover, in order to select the VBF 
contributions, the two tagging jets are required to reside in opposite detector hemispheres 
and to be well separated in rapidity, 

Vji > > r/j2, Arjjj = rjj^ - r]j^ > Ar]jj^in. (5.2) 

Varying the value of Ar/jjmin, we study the fraction of the VBF contributions to the cross 
section with the exact matrix elements. The analyses are done at the parton level with 
tree-level matrix elements. The exact matrix elements for the Higgs boson productions 
are calculated by HELAS subroutines [38], generated by the HEFT (Higgs effective field the- 
ory) model in MadGraph/MadEvent (MG/ME) v4 [48]. For the massive-graviton produc- 
tions, the relevant HELAS subroutines for massive spin-2 particles and its interactions based 
on the effective Lagrangian of eq. ( fl.30| ) have also been implemented into MG/ME [46]. 
Numerical integrations are done with the help of the Monte Carlo integration program 
BASES [49]. Throughout our numerical study, we employ the CTEQ6L1 parton distribu- 
tion functions [50] with the factorization scale chosen as the geometric mean of the jet 
transverse momenta /jf = y/Pf~PT~, and fix the QCD coupling at = a^^{mz) = 0.13. 
Unless specified, we set the heavy particle mass at M = 600 GeV and the scale of the RS 
model at A = 4 TeV, which corresponds to the current lower bound for the mass of the first 
KK mode of massive gravitons [17, 24]. Note that we use the constant value for the Higgs 
effective couplings in spite of Mh,a > 2mt, because the energy dependence of the effective 
couplings are almost canceled out when the ratio of the cross sections is considered, and 
because the azimuthal angle distributions are insensitive to the form factor effects [8]. 

In table 0, we show the ratio of the VBF contribution (ctvbf) to the cross section with 
the exact matrix elements (cxexact) for the nine subprocesses at the LHC, after imposing the 



inclusive cuts (5T) and the VBF selection cuts (|5.2|) with Arjjj^i^ = 3, 4 and 5. For all the 
subprocesses, as the rapidity separation increases, the VBF contributions tend to dominate 
the exact matrix elements and the ratios approach unity. In the first row, the CP-even 
and -odd Higgs production processes in qq collisions, where we consider ud collisions for 
simplicity, have only the VBF diagrams, and hence cvbf /inexact = 1.^ Although the cross 
sections in qg and gg collisions have non-VBF diagrams, their contributions are rather small 



^Even for the collisions of identical-flavor quarks, the s- or ii-channel contribution is negligible when the 
VBF cuts are applied [11]. In addition, the interferences between the electroweak and QCD contributions 
are very small [12, 51, 52]. 
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C^VEF exact 


Arijj > 3 Arijj > 4 A??,,- > 5 


qq ^ qqH/A/G 
qg ^ qgH/A/G 
gg ^ ggH/A/G 


1.00/1.00/1.02 1.00/1.00/1.02 1.00/1.00/1.02 
1.04/1.04/1.07 1.03/1.03/1.06 1.02/1.02/1.04 
1.05/1.05/1.09 1.04/1.04/1.07 1.02/1.02/1.05 



Table 8: The same as table @, but imposing the additional ptj cut of eq. (|5.3| ). 



after the VBF cuts. On the other hand, for the massive-graviton production, the gravitons 
are emitted also from the quark lines, as shown in figs. ^ (a) and (b), and the ratio deviates 
significantly from unity especially in qq collisions. Furthermore, their contributions do not 
diminish swiftly when the rapidity separation cut is increased from 3 to 5. We note here 
that the electroweak contributions to the qq — > qqG process represent a small correction, 
which is below 1%, even when the VBF cuts are imposed [27]. 

In case of the massive-graviton production, the non-VBF contributions are significant 
even after the VBF selection cuts with Arjjj^i^ = 5. We therefore examine the impact of 
an additional cut on the transverse momenta of the tagging jets, 



PT, < 100 GeV. 



(5.3) 



Table ^ shows the same ratio of the VBF contribution to the exact cross section when this 
additional cut is imposed. We find that the above pTj cut works effectively to suppress 
contributions from the non-VBF diagrams, especially the diagrams which emit the graviton 
from the quark lines. This is because the quarks that emit a graviton tends to have high 
transverse momenta. 

Summing up, from table ^, the large rapidity separation may guarantee the validity of 
our VBF analyses not only for the WBF Higgs productions but also for the GF processes. 
Moreover, from table H], the pxj slicing cut. 



20 GeV < pT. < 100 GeV, 



(5.4) 



from eqs. ( |5.1| ) and (5.3) is effective in selecting the GF contribution to the graviton pro- 
duction processes. It should be noticed that stringent cuts increase the VBF contributions 
but reduce the primary event number. 

As a reference, we present the total cross sections with the exact matrix elements for 
the massive-graviton productions in table ^ where the inclusive cuts ( |5.1| ) and the VBF 



O-exact [Pb] 


20 GeV < pr, 
Ar?j, >3 Ar/,-, >4 Ar/,-, > 5 


20 GeV < pT, < 100 GeV 
Ar/,-, >3 A77,-, >4 Ar/,-, > 5 


qq qqG 

qg qgG 
99 ggG 


2.0 1.6 1.1 
13.2 8.8 4.7 
15.9 8.0 3.2 


0.7 0.6 0.5 
7.2 5.5 3.3 
11.7 6.2 2.7 



Table 9: The total cross sections with the exact matrix elements for each subprocess of the massive- 
graviton plus dijet events at the LHC, after imposing the inclusive cuts ( |5.lD and the VBF cuts (5.2) 
without and with the additional pt, cut of eq. (5.3), for Mq = 600 GeV and A = 4 TeV. 
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cuts ( |5.2| ) without and with the shcing cut (|5.4| ) are imposed. For the cross section 
of the Higgs boson productions, see fig. 4 in ref. [2]. In the following analyses, we take 
Arjjjaiin = 4 in the VBF cuts ( |5.2| ) for the Higgs and graviton productions, and further 
apply the pT slicing cut (5.4) for the graviton productions. 



5.2 Correlations in the production with two associated jets 

We are now able to discuss the angular correlations of the two accompanying jets in the X 
productions, using our analytical VBF amplitudes. 



The production and decay density matrices in eq. (2.17) have all information on the 



angular correlations between the jets in the productions and the decays of the heavy particle 



X. In this section, we consider the n = 2 case in ( 2.17 ) to investigate the angular correlation 
between the two jets in the Xjj productions, while, to simplify the decay part, we consider 
the n' = 2 case and fix the X polarization along the 2-body decay axis (z'-axis). A' = 
a'l — a'2 = ±2, ±1, or 0. We note that, in practice, we can project out a'l and erg by 
properly weighting the final states of a'^ and decays in the X — > a'^a^ decays; e.g. for 
the X W~^W~ decays we can project out all the five cases. A' = ±2, ±1, and 0, while 
the X ^ 77 or gg decays give only the sum of A' = +2 and —2. Moreover, the X t^t^ 
decay process can project out A' = +1 and —1 cases, while for the e or ^ case we cannot 
distinguish these two. In the following we take A' = A' and suppress the decay density 
matrix in eq. ( 2.17| ). 



For the VBF processes (|2.1| ), the production density matrix is given in eq. ( 2.19a| ) in 



terms of the production tensor V^^^ of eq. ( p. 20a ). The tensor has 81 independent jet 



angular distributions in terms of the polar (^1,2) and azimuthal ((/>i,2) angles of the two 
tagging jets. When we isolate the azimuthal angle dependence in eq. ( |2.17 ), there are 25 
distributions (including one constant piece) as 

V^^J'^^diydl ^, = + {2'Ste [F^ cos </>i + F3 cos (/>2 + F4 cos 2(/>i + F5 cos 2(/>2 

+ F^ cos((/)i ± 02) + F^ cos(20i ± (1)2) + F^ cos(0i ± 2^2) 
+ ^9^ cos 2(01 ±02)] +(3^e — > 9m, cos —>■ sin)}. (5-5) 

Here, and in the following, summation over repeated indices (Ai, A2, Ai, A2) = ±,0 is im- 
plied, and a shorthand notation such as F^ cos(0i±02) for Fq' cos(0i+02)+-^^6" cos(0i— ^2) 
is used. The coefficients -F/^^ are the functions of the kinematical variables except the az- 
imuthal angles 0i,2- For the productions of spin- full heavy particles, they also depend 
on the decay angle G which comes from the product of two d functions. For the spin-0 
particle case, only the five terms in eq. ( |5.5D survive due to the helicity selection Ai = A2, 



which will be discussed in the next subsection 5.2.1. All the sine terms vanish when CP 
is conserved and when the absorptive part of the amplitudes are neglected, e.g., in the 
tree-level approximation. It should be noted that the azimuthal angle variables, 0i and 
02) are individually defined in the VBF frame by the scattering plane of the subprocess, 
V*V^ ^ X ^ a\a'2] see also fig. |. 

Because the phases of the quark and gluon current amplitudes are the same for each 
helicity combination (see tables |l|(top) and ^(top)), and because the phase of the product 
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of the two currents for ai = 0"3 and (T2 = (T4 is (see eq. ([4.6D ) 

Jlt^,=^j2t<r,=a, OC e-(^l^-^^^2) , (5.6) 

the coefficients F^^g for (7(7, qg and (7(jr collisions are expressed in terms of the production 

tensors Vt^t^ and two d functions as 

A1A2 

-^3 = ^Ai^+'^Ai,A'^Ai-l,A' + ^AiO "^Ai+l.A'^^Ai.A'! 
-^4 = 'P-A2'^l'-A2,A''^-1-A2,A'' 
-^5 = '^Ai+'^Ai+l,A''^Ai-l,AM 

= 'Po±di,y^-l/l,X' + ^00''' ^2/0, A' ^0, A' + ^-±^0,A'^-2/0,A' + ^-0^1/-1,A'^-1,A') 

= 'P^±^l,A'^-2/0,A' + ^^cr'^2/0,A'^-l,A'> 
^8 = ^Oi'^VCA'^-l/l.A' + ^-±^l/-l,A''^-2/0,A'' 

-^9" = ^-±^2/0,A'^-2/0,A'- (^•'^) 

This relations can be also applied to the processes with the helicity-flip currents, namely the 
cases for cJi = — CJ3 and/or CJ2 = — 1T4, which appear only in the gluon currents, although 
they have the different phase from those with the conserved currents in eq. ( |5.6D . The 
process that one of the currents is helicity-flip (the case in ai = —cr^ or (72 = — CT4) leads 
nontrivial azimuthal distributions from F2 through F^, while the case that the both currents 
are helicity-flip gives rise to only the constant piece, Fi. This is because the helicity-flip 
gluon splitting amplitudes emit an off-shell gluon of definite helicity; see table ^. On the 
other hand, the helicity-conserved processes (ui = (T3 and (T2 = cr^) can contribute to all 
the terms, as can be seen from table |l] for the quark currents and table § for the gluon 
currents. We note that the magnitude of the correlations is determined by the relative 
ratio to the constant term Fi. 



It should be stressed here that in eq. (5.7) the ©-dependent azimuthal angle correlations 



are manifestly expressed by quantum interference among different helicity states of the 
intermediate vector-bosons. We also notice that the above formulae can be applied to any 
spin-J particle productions through the VBF processes, although massive spin-0 and -2 
particles are considered in this article. 

5.2.1 Higgs boson productions 

For the scalar particle productions, only the three off-shell VBF amplitudes, Aix++, ■Mxoo 

and Mx 1 in which the colliding vector-bosons have the same helicities (Ai = A2), 

can contribute to the production amplitude (2.7), and there is no Q dependence, namely 



d'l = 1 in eq. ( p.5| ) (and in eq. (p. 7])). Therefore, the azimuthal angle correlation (| 

is reduced to 

p^i^2 =F^ + {2?fie[FQ cos A(?ii2 + Fq cos2A(/)i2] + (5fte ^ 9m, cos ^ sin)} (5.8) 
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with A(/>i2 = (j)! — 4>2, which is the azimuthal angle separation of the two tagging jets. 
As mentioned above, the azimuthal dependence is manifestly expressed by the quantum 
interference terms among different helicity states of the intermediate vector-bosons; the 
Fq (= Vqq^ + V^_) term is induced through the interference between the production 
with the longitudinal and transverse polarization states of the vector bosons, while the 
Fg (= V^^) term is caused by the interference between the two transverse polarization 
states. 

Now we can observe clearly the origin of the nontrivial azimuthal angle correlations 
for the CP-even and CP-odd Higgs bosons {H and A), predicted in refs. [1~4]. As we 



mentioned in section [4.3| , in the case of Qi,Q2 ^ M, where the VBF contributions are 
dominant, the Higgs bosons with a g^'^-type coupling are produced mostly through the 
longitudinally-polarized vector-bosons. Therefore, the Mhoo amplitude dominates the 
total amplitudes, and hence there is little interference terms in eq. This is why the 

WBF processes give the flat azimuthal angle correlation, 

d^H(WBF) ~ ^1- (5.9) 

For the loop-induced GF Higgs boson couplings, on the other hand, they are mainly 
produced by the transversely-polarized vector-bosons, namely the M.h/a++ -^H/A 



amplitudes have the dominant contribution. Therefore, by the relations in (4.32) and 
( |4.33 ), the azimuthal distributions are 

dan/AiGF) ~ [Pi ± 2 \Fg \ cos 2A^i2\ , (5.10) 

where the +/— sign is for CP-even/odd Higgs bosons. One can clearly see that the 
azimuthal distribution is strongly suppressed (enhanced) around A4>i2 = 7r/2 for the GF 
CP-even (-odd) Higgs boson productions. We note that the Pg~ term in eq. ( |5.8| ) is exactly 
zero not only for the CP-odd Higgs boson production but also for the GF CP-even Higgs 
boson production since it measures the P-odd amplitude. 

From the relations in and by using the explicit forms of the amplitudes in ta- 

bles ll(top), §(top) and ^ the coefficient functions Pi and Pg~ in eq. ( 5.101) for each sub- 
process, qq, qg and gg scatterings, are given by 

Fi[qq] = (l + cos2 0i)(l + cos2 02), (5.11a) 

= ^--(l + cos2 0i)(l + 3cos2 02)^ or (1^2), (5.11b) 
sm 82 

= -^^^-^(l + 3cos2 0l)2(l + 3cos2 02)^ (5.11c) 
sm 01 sm 02 

and 

Fg[qq/qg/gg] = ±^{l - cos^ ei){l - cos^ 02). (5.12) 

Here we take the Qi, Q2 M limit for the VV X amplitudes, where the only surviving 
amplitudes are = (see table ^), with the common overall factor 

(±) ^ atgjj/AggM^QlQl 

' 2cOs2eiCOs2^2 ' ' ^ ^ ^ 
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Note that we suppress the color factors, which are relevant to the ratio of the qq, qg and 
gg contributions in realistic simulations. It is remarkable that the interference term Fg~, 
which receives the contribution only from the helicity-conserved amphtudes (cJi = (T3 and 
CJ2 = C74), is same for all qq, qg and gg collision processes. Meanwhile Fi has the different 
contributions from the quark currents and the gluon currents, which includes the helicity- 
flip contributions (cJi = —as and/or a2 = — o"4), and can be larger as the process involves 
the gluon currents. These indicate that the gluon currents reduce the interference effect. It 
may be worth presenting the above functions in terms of the 2:1^2 variables in the collinear 
limit (/3i,2 1) in eq. ( ^Isl) : 



and 



with 



Fi[qq] = 4(2 - 2zi + zf)(2 - 2z2 + 4)> (5-14a) 

Fi[qg] = —^-A2-2zi+zl){l-Z2 + zlf or (1^2), (5.14b) 
(1 - Z2) 

Fi[gg] = J, r(l -zi+ zlf{l -Z2 + zlf, (5.14c) 

(1 - zi)y^ - 



F^[qq/qg/gg] = ±8(1 - z^){l - Z2), (5.15) 



2z^ Z2 

To examine the validity of the above analytic parton-level expectations, we plot in 
fig. ^ the normalized azimuthal correlations A(f)i2 (mod 27r) between the two tagging jets 



in the Higgs + 2-jet productions at the LHC, where the selection cuts (|5.l|) and (|5.2D with 
Aryjjmin = 4 are applied. The distributions for each subprocess with the full diagrams and 
those with the VBF diagrams only are shown by solid and dashed lines, respectively. The 
VBF contributions can reproduce the distributions with the exact matrix elements very 
well not only for the WBF processes but also for the GF processes. As mentioned in the 
introduction, these azimuthal angle correlations predicted in the leading order may survive 
even after higher-order corrections are applied [13, 15]. 

It may be worth pointing out here that, even though our definition of the azimuthal 
angles of the jets, which are measured along the vector-boson colliding axis, is different from 
the usual definition along the beam axis in the laboratory frame, the A(/>i2 distributions 
are almost same in the two frames due to the characteristic VBF kinematics. The VBF 
amplitudes are dominant in the collinear limit ( [4.181 ), where the vector-boson colliding axis 
(i.e. the z-axis in fig. ^) is identical with the beam axis, and hence the A(pi2 distributions 
are not so much distorted. In fact, our A(j)i2 distributions for the Higgs boson productions 
semi-quantitatively confirm those in the previous works [1-4]. 

Before turning to the spin-2 case, there are a few remarks related to the previous 
studies on the azimuthal correlations in the Hjj events, (i) The XVV coupling form 



factors are factorized as in eq. ( 5.13| ), and therefore, the A(j)i2 distribution is insensitive to 
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(a) qq^qqX (b) qg^qgX (c) gg^ggX 




A4>i2 = 4>i - 4>2 



Figure 5: Normalized azimuthal correlations A0i2 (mod 27r) between the two tagging jets in 
the pp — > jjX process at the LHC, where the selection cuts (5^) and ( ^.21) w ith Ar]jjaun = 4 
are imposed. For the massive-graviton productions, the additional pr^ cut (^j) is also imposed. 
The distributions for each subprocess with the full diagrams (solid lines) and with the only VBF 
diagrams (dashed lines) are shown. 



their effects [8]. (ii) As the ratio of the Higgs boson mass M to the partonic CM energy 
\/I decreases, the interference effect grows [8]. In the collinear limit ( 4.18 ), we obtain 

s 



-4 cos 9i cos 6*2 



cos 02 



(5.17) 



(l + cos6li)(l 

Therefore, as M/^/I becomes smaller, i.e. z\Z2 0, cos^i and cos 02 approach zero, and 
the ratio of Fg~ to F\ in eq. ( ^.10 ) grows; see eqs. ( ^.11 ) and ( ^.121) . (iii) Although we have 
considered the three types of tensor structures separately for the Higgs coupling to vector 
bosons in this paper, it is easy to extend our analyses to a mixed CP scenario [3,4, 10]; 
for instance, the additional phases, which come from the CP-mixed XVV coupling, can 
give rise to the sine terms in eq. ( |5.8D , and explain the shift of the dip positions in fig. 8 of 
ref. [4]. 

5.2.2 Massive graviton productions 

Here, we discuss the case for the spin-2 particle productions, which is more involved than 
the scalar case, because all the nine amplitudes generically contribute to the total amplitude 



in eq. (2/7), which can lead all the 25 azimuthal distributions in eq. (|5.5| ). Moreover, the 
graviton polarization along the momentum direction of the decay products (A') depends 
on the decay angle 0. 



In section 5J. we demonstrated that the QCD VBF amplitudes can have significant 
contribution to the G + 2-jet events by imposing the VBF cuts and the pr, slicing cut. 
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In this case, the two off-shell VBF amplitudes, A4g+- ^-nd AicZ^, are dominant, as 
mentioned in section |4.3| . Therefore, only the (oc V^^) term in eq. ( |5.5| ) dominantly 
gives the nontrivial azimuthal correlation, 



d&G ~ [Fi + 2 F9+ cos 2$ 



12 



(5.18) 



with ^12 = 4>i + 4>2- It should be emphasized here that <&i2 is not the azimuthal separation 
^012 (= 01 — (P2), but the sum of the azimuthal angles of the two jets, <j)i and (j)2- 

From ( |5.18 ), one can immediately conclude that the A(j)i2 distributions for the massive- 
graviton productions are flat. In fig. ||, the A(pi2 correlations for the KK graviton produc- 



tions are also plotted with thick lines, where the pT slicing cut (5^) has been imposed 



as well as the inclusive cuts (5J.) and the VBF cuts ( ^.2[ ).° The contributions from the 
A = ±1 and states, which can give rise to the cosA(/)i2 and cos2A0i2 dependence in 
eq. ( ^.SD , are invisibly small; they are smaller than the A = ±2 by two and three orders of 
magnitude, respectively. The flat A(f)i2 distribution for the massive-graviton productions 
is distinct from that for the SM Higgs boson productions, which is expected to have a dip 
around A(j)i2 = n/2 due to the GF contributions. 

Now let us see the explicit forms of the functions Fi and in eq. ( |5.18D . From 
eq. (5/7) and by using the explicit forms of the amplitudes in tables ||(top), §(top) and |5|, 
one finds 



F^iqq/qg/gg] = F^[qq/qg/gg] x \{{dl^^AQ)Y + {d\AQ)f}, 



F 



-H 



qq/qg/gg] x dl2v{@)d\y{e), 



(5.19) 
(5.20) 



where F^ and Fg ^ are the same as in eqs. ( 5.11 ) and ( ^.12 ), respectively, for the Higgs 
boson productions, and the common overall factor is 



COS^ 61 COS^ 6*2 



(5.21) 



Here, from table |5|, we take A^c^^ = — in the collinear limit. By explicit forms of 
J = 2 d functions (see table ^), the functions in eqs. ( p. 19 ) and (5.20) are 



'^(1 + 6cos2e + cos^e) for A' = ±2, 



F^[qq/qg/gg] = f^ [qq/qg/gg] x < 



i(i -cos^e) 



I sin^ G 



for A' 
for A' 



±1, 
0, 



(5.22) 



and 



F^^[qq/qg/gg] = Fg "[qq/qg/gg] x < 



^ sin^ G for A' = ±2, 
-isin^G forA' = ±l, 
I sin^ G for A' = 0. 



(5.23) 



*See also fig. 2 in ref. [46], where Mq ~ 1 TeV and the different selection cuts are applied. 
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Figure 6: Azimuthal correlations <i>i2 (mod 2tt) between the two tagging jets for the gg-scattering 
subprocess in Gjj productions at the LHC, where the same selection cuts in fig. |^ are imposed. The 
distributions for each final G polarization state, A' — ±2 (a), ±1 (b), and (c), are shown at the 
decay angle Q — n/2 (dashed), 7r/3 (dashed-dotted), and n/6 (dotted), where they are normahzed 
by a{pp — !■ jjG) B{G a-'ia'2)- The distributions after integrating out Q are also shown by thick 
solid lines. 

The distributions strongly depend on the decay angle and the final polarization A'. The 
dependence is the same for A' = +2 and —2, and also for A' = +1 and —1. Since 
J = resonances do not have such dependence, the J = 2 and A' = state can be 
distinguished from the J = A' = state in principle. The coefficient function for all 
the final polarization states is proportional to sin^ G, and hence at = and vr, where the 
decay axis (z'-axis) is coincide with the initial polarization axis of the gravitons (z-axis), 
the azimuthal correlation is absent. This is because each colliding vector-boson has the 
definite helicity. Meanwhile, the correlation becomes larger with the larger decay angle, and 
reaches the maximum at = 7r/2, where the A = +2 and —2 states are mixed maximally. 
It should be noted that the sign of the Fi^ term for A' = ±1 is different for the other states, 
which gives rise to the distinctive correlation. 

After integrating out the decay angle 0, we obtain 

Ff'iqq/qg/gg] = F^[qq/qg/gg\ X ^{1, 1, 1} for A' = {±2, ±1,0}, (5.24) 

F^^iqq/qg/gg] = F^"[qq/qg/gg] x ^{i, -4, 6} for A' = {±2, ±1,0}. (5.25) 

The above results show that the magnitude of the azimuthal correlation depends on the 
final polarization A'. The interference effect for the A' = case is largest among the possible 
five polarization states, although the decay branching ratio of KK gravitons into the A' = 
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state, such as a pair of longitudinal weak-bosons and a ti pair with same helicities, is 
less than 1% for a whole mass range of gravitons. The correlation for A' = ±1, which is 
realized in the decays into a fermion pair, is four times larger than that for A' = ±2, which 
is projected out in the G — > VV decays. It is worth noting that the above results depend 
only on A', not on the decay mode. This universality of the angular correlation can be an 
experimental signal of the X spin measurement. 

To examine the above analytic expectations, we demonstrate in fig. ^ the azimuthal 
correlations $12 (mod 2-ir) between the two tagging jets for the g(;-collision subprocess in 
Gjj events at the LHC, where the full diagrams are taken into account and the same 
selection cuts in fig. ^ are imposed, i.e., the selection cuts (5J), ( |5.2| ) and ( |5.3D . The 



distributions for each final polarization state. A' = ±2, ±1, 0, are shown at the decay angle 
= tt/2, 7r/3, and 7r/6, where they are normalized by a{pp jjG)B{G — > a'ia'2). The 
distributions after integrating out G are also shown by thick solid lines, which is normalized 
to unity. The above analytical results can describe the simulations very well. For the qg- 
and ^(^-scattering cases, all the qualitative behaviors are same, but the interference effects 
diminish as in the A(/>i2 correlations for the Higgs boson case in fig. |5[ It should be stressed 
here that the observation of the d function behavior, or the @ dependence, is a measurement 
of the X spin, which can be strengthened by the azimuthal correlation between the tagging 
jets. 

5.3 Correlations in the decay into a vector-boson pair 

As mentioned in section |2|, the processes of a heavy-particle decay into a vector-boson pair 
which subsequently decay into M, qq, or gg are closely related to the VBF production 
processes. Here we discuss the X decay correlations between the two decay planes of the 
vector bosons by using the explicit helicity amplitudes, as in the production process. 

To simplify the production part, we consider s-channel X productions in gg fusion or 
qq annihilation and its subsequent decays into four-body final states, namely n = and 
n' = 4 in (2.17), so that the initial polarization of X along the z-axis can be fixed, A = ±2, 



lb 1, or 0. In the following we take A = A and suppress the production density matrix in 
eq. 

For the X decay processes ( 2.10| ), the decay density matrix is given in eq. ( 2.19b| ) in 



terms of the decay tensor T^^J^J of eq. (2.20b). Similar to the azimuthal angle distributions 



for the production in eq. ( |5.5D , those for the decays into a vector-boson pair in eq. ( 2.17| ) 
are generally expressed by [53] 

diyd-l -yV^}^f =F[ + {2^e [F^ cos (/)[ + F3 cos 02 + cos 2^; + F5 cos 2(j)'2 

+ F^^ cos{(p[ ± (p'^) + F!^^ cos{2(l)[ ± (p'^) + F^^ cos(</>; ± 2(^'2) 
+ Fg^ cos 2{<j3[ ± (p^)] + (5Re ^ 9m, cos ^ sin)}. (5.26) 

Here, and in the following, summation over repeated indices (A'j^, A2, A'^, A2) = ±,0 is im- 
plied, and 4>^ is the azimuthal angle between the decay plane of the vector boson and the 
X production plane {gg/qq ^ X — > VV) in the partonic CM frame; see also fig. |3[ 
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Because the azimuthal angle dependences of the quark and gluon current amphtudes 
are the same for each hehcity combination (see tables |l](bottom) and ^(bottom)), and 
because the phase of the product of the two currents for a'l = — CJ3 and a'2 = —cr'^ is (see 
eqs. (pip and ( p^ )) 



,J6„'- oce^(^'i<^'i~^^^<^^), (5.27) 
the coefficients F'-^^}^ for (//)(//), {ff){gg) and {gg){gg) decays are given in terms of two 



A' A' 

d functions and the decay tensors T>^}^J as 

7^ rlJ Tl^'lK 



H = ^^A -A^^A,1-AJ,^+a| + ^A,-1-A'2^A -AJj^OA'a^ 
-^3 = ^A,A'i-l^A,A'j^A'jO^ + ^A,a;^A,a;+1^a;-' 

— "a -1-A^"A,1-A^^+A^' 
-^5 = ^A,A'i-1^A,A'i+1^a(-' 
-^6^^ = ^A,-l/l^A,l^+0 + ^A,0^A,2/0^+T + '^A -2/O^A,0^od^ + '^A -1^A,1/-1^0=F ' 
= ^A- 2/0^11^+0" + ^A,-l^A,2/0^+T' 
= ^A,-l/l^A,2/0^+T + ^A,-2/0'^A,l/-1^0=F^' 
-^9"^ = ^A -2/0^A,2/0^+T- {^.2^) 

We note that this relations can be also applied to the processes with the same-helicity 
currents, namely the cases for cr^ = a'^ and/or = C4, which appear only in the g* — > gg 



currents; the same arguments in the production part can be applied (see below eq. (5.7)) 



We should note that, similar to eq. (5^) for the VBF production processes, the 0-dependent 
azimuthal angle correlations are manifestly expressed by quantum interference among dif- 
ferent helicity states of the intermediate vector-bosons in eq. ( ^.28 ). 



5.3.1 Higgs boson decays 

The decay amplitudes ( 2.15| ) for scalar particles are the coherent sum of the three ampli- 



tudes in which the decaying vector-bosons have the same helicities (A'^ = A2), and there is 
no G dependence, i.e. d^^?(0) = 1 in eq. (|53^ ) (and in eq. ( ^.28D ). Therefore, similar to 



the production case, the azimuthal angle correlation ( ^.26 ) for the J = particle decays is 
reduced to 

pjijf =F[ + {2^e[FQ- cos A(/.;2 + coslA^'^^] + (3f?e ^ 9m, cos ^ sin)}, (5.29) 



where A(/>'^2 = 'P'l ~ ^2 is the angle between the two decay planes of the vector bosons. 

The decay angular correlations have been extensively studied for the Higgs boson 
decays into the weak bosons, B. Z^*^ Z , W^*'^W [29-33]. It may be worth pointing 
that the Fq~ and F^~ terms in eq. ( ^.29 ) can be observable for the Higgs boson mass 



-33- 



Mh < 400 GeV since the decays into a pair of transverse weak-bosons are not negligible 



compared to the decays into longitudinal ones, as shown in (4.34). Fq~ is very small for the 
H — > ZZ decay, while it is large for the H WW decay due to the P-odd nature [30, 34]. 
One can see detailed simulations for the LHC in refs. [31-33]. It must be recalled here that 
there is no azimuthal correlation between the two jets in the WBF Higgs productions as 
discussed in eq. (5.9). 



Here we discuss the correlations for the Higgs boson decays into a virtual-gluon pair, 
H/A — > g*g* — > {qq){qq), {qq){gg), or {gg){gg), in some detail. In the on-shell limit, or 
Q'i 2 0, where the above decay processes via a virtual-gluon pair are dominant, only the 
transverse amplitudes {M.'jj^^±±) contribute to the decay amplitudes both for the CP- 
even/odd Higgs boson decays. Moreover, due to the P-odd nature, the Pg~ term vanishes 
for the QCD processes as in the VBF processes. Therefore, only the Fg~ term in eq. ( 5.2S| ) 
is relevant in this case. 

From the relation in (|5.2g| ) and by using the explicit forms of the amplitudes in ta- 
bles |l|(bottom) , |2|(bottom) and ^, the coefficients F[ and Fq~ in eq. ( 5.2g| ) for each decay 



process, {qq)iqq), {qq){gg) and {gg){gg), are expressed (except for the color factors) as 

H[{qq){qq)] = (l + cos2 0;)(l + cos2 0^), (5.30a) 

F[[{qq){gg)] = -^{l+cos^e[){3 + cos^e'2f or (1^2), (5.30b) 
sm 6*2 

F[[{gg)igg)] = -^-^-^(3 + cos^ 0;)2(3 + cos^ 0^)2, (5.30c) 
sm (J\ sm 

and 

H~[{qq){QQ)/{9g){99)] = ±^(1 - cos^O[){l - cos^9',), (5.31a) 

F^-[{qq){gg)] = Tl{l-cos^0[){l-cos^9',), (5.31b) 

where the upper (lower) sign corresponds to the CP-even (-odd) Higgs boson case. Here 
we take the ^ limit for the X VV amplitudes, where the only surviving 

amplitudes are |A^^^^±±| = ^^^1 with the common overall factor 

Fl^^^ = \9\iA,/sM'Q?Q? F'X^\ (5.32) 

The constant term F'^ has the different contributions from the qq decay and the gg decay and 
can be larger as the process involves the g* gg splitting. As in the production processes 
with the helicity-flip currents, the same-helicity currents {a'^ = cr'^ and/or = cr'^), which 
appear only in the gluon currents, give rise only the constant piece, F[. The interference 
term Pg", which receives the contribution only from the opposite- helicity currents {a[ = 
— CT3 and a'2 = —'74), is same except for sign in all the decay processes. This is because the 
relative sign between the different helicity states for the outgoing fermion and gluon current 
amplitudes is different; see tables |l](bottom) and ^(bottom), and because the relative phases 
between Jl^/ and j!^, for the quark currents ([4.11| ) are different from those for 



the gluon currents (4.15) 
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Aside from the coupling constant, the double (qq) decays via a pair of gluons are the 
same as the decays via a virtual-photon pair, H/A — > 7*7* — > {U/qq)[U/qq)^ and we expect 
the strong correlations as in the GF qq qqH/A processes. On the other hand, the decay 
correlations for {qq){gg) and {gg){gg) are much smaller than those for qg — > qgH/A and 
gg ggH/A since the reduction factor from the gluon currents for the decays, (3+cos^ O'-)"^ 
in ( |5.30| ), can be much larger than that for the productions, (1 + 3cos^ in ( |5.11| ). This 
is well known as the knowledge of the QCD parton branching; a quark-antiquark pair from 
the gluon, g* — > qq, has the strong correlation between the gluon polarization and the 
decay plane, while for g* gg the correlation is weak [54]. It may be a challenging task 
to observe the azimuthal correlations between the QCD jets in the Higgs boson decays at 
the LHC. The quantitative study will be reported elsewhere. 



5.3.2 Massive graviton decays 

For the spin-2 particle decays, unlike the scalar particle decays, the nine amplitudes should 
be coherently summed in eq. ( 2.15| ) and may give rise to all the 25 azimuthal distributions 



in eq. §2§ 



In this article, we study the process of the graviton decays into a weak-boson pair at 
the LHC, 

pp^G^ WW/ZZ ^ (/i/3)(/2/4). (5.33) 

This includes final four-charged-lepton signals, which may give a clean signal and allow 
a complete kinematical reconstruction at the LHC. The s-channel G production has two 
possible sources, gg fusion and qq annihilation, even though the gg contribution dominates 
the cross section for the graviton mass up to 3.4 TeV [19]. In the parton CM frame, the 
polarization of the produced gravitons is fixed along the beam axis (z-axis) at A = +2 or 
—2 in gluon fusion, while at A = +1 or —1 in qq annihilation. The two different production 
modes lead totally different angular distributions, as we will see below. 



As mentioned in section O, in the heavy graviton mass limit (/?' = y 1 — 4m'y/M'^ —>■ 
1), only three G VV amplitudes survive for the decays into on-shell weak-bosons; see 
table |6|. In this limit, therefore, the azimuthal distributions in eq. ( |5.26| ) can be reduced to 

dl^^,dl -yV^}^f =F[ + 2 cos $'12 + 2 i^9+ cos 2^'^^ (5.34) 

with ^'^2 = 't^'i + 02- It should be stressed here that ^'^2 is not the angle between two decay 
planes A(/)']^2 (= ~ ^2)' sum of (j)'i and (/>2, which are measured separately from 

the graviton production plane. 

From eq. ( |5.28D and by using the explicit forms of the helicity amplitudes in ta- 
bles ll(bottom) and |^ and d functions in table |l^, the complete angular distributions of 
eq. ( |5.34| ) are presented for A = ±2 and ±1 in appendix 0. After integrating out 6'^ and 
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^2 in ( P-l|) and ( C.2 ), the 0-dependent coefficients in eq. ( ^.34|) are 





= 3 + 10 cos 




97r^ 








= isin^e, 



^G + 3cos^G, (5.35a) 
1-cos^G), (5.35b) 

(5.35c) 



for A = ±2 (the gg initial state), and 



F{ = 8 + 4cos2G - 12cos'^e, (5.36a) 

F'+ = -KiK2 sin^ G cos^ G, (5.36b) 

16 

F9+ = -sin^Q, (5.36c) 
for A = ±1 (the qq initial state), where the common over all factor is 

Fl^^^ = ^9hvv9h.h9h.uM'Q'M' Fl(^^ (5.37) 
with gy^f.f.^.^ = i^g^^^^'-+'^Y + {g+^^^^^^Y , and the combination of the V f f couplings Kj is 

The V f f couplings are given in ( |4.3| ). The G distributions for the $'^2"i'^dependent term, 
F[, are totally different between the two initial states due to the angular momentum con- 
servation; a pair of vector bosons from the A = ±2 gravitons tend to decay to the forward 
and backward directions, while the A = ±1 gravitons are not allowed to decay into a 
vector-boson pair at G = and vr (in the /?' = 1 limit), where the momentum direction of 
the decaying vector-bosons (z'-axis) is coincide with the polarization axis of the gravitons 
(z-axis). On the other hand, all the azimuthal correlations, Fg"*" and -Fg"*", disappear at 
G = and vr since each decaying vector-boson has the definite helicity. The coefficient Fg"*" 
takes the maximum at Q = 7r/2 (vr/4) for A = ±2 (±1), while the Fg"*" becomes larger with 
the larger decay angle and reaches the maximum at G = 7r/2 both for A = ±2 and ±1. 
The sign difference of the $'^2"dependent terms between A = ±2 and ±1 gives rise to the 
distinctive correlations. 

Let us estimate the asymmetries Ai = 2Fl^/F[ {i = 6,9) from ( 5.35| ) and ( 5.36| ) for 
the G — > ZZ M process. After the integration for G, one finds 

Aq = KiK2 0.007, Ag = 0.038 for A = ±2, (5.39) 

832 26 

97r2 2 

Aq = -KiK2 0.005, ^9 = 0.154 for A = ±1. (5.40) 

416 13 

The asymmetry Aq of the cos<I>']^2 term is tiny both for A = ±2 and ±1, less than 1%, due 
to the smallness of the parity violation for the Z decays, similar to the H — > ZZ M 
process. On the other hand, Aq of the cos 2$'^2 term reaches around 4% for A = ±2 and 
15% for A = ±1. 
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Figure 7: Normalized azimutlial correlations (mod 27r) between the decay planes of the vector 
bosons in the G ZZ M decays for the initial helicity A — ±2 (a) and ±1 (b), where the decay 
angle 8 is fixed at 7r/2 (dashed), 7r/3 (dashed-dotted) , and 7r/6 (dotted). The distributions after 
integrating out 6 are also shown by thick solid lines. 



We simulate the G — > ZZ — > 4£ process to examine our analytic expectations, where we 
take M = 600 GeV, mz = 91.2 GeV and Tz = 2.5 GeV, that is, all the nine amplitudes in 
(2.15) are taken into account without any approximation. Figure ^ shows the normalized 
azimuthal correlations ^'12 (mod 2ti) between the decay planes of the vector bosons for 
A = ±2 (a) and ±1 (b). The decay angle is fixed at 7r/2, vr/S, and 7r/6, shown by dashed, 
dashed-dotted, and dotted lines, respectively. The distributions after integrating out G are 
also shown by thick solid lines, which are normalized to unity. Our analytic approximation 
can explain the results not only qualitatively but also quantitatively. It should be stressed 
again that the observation of the dependence of the azimuthal correlations can be a 
measurement of the X spin. 

At the LHC, as mentioned before, gg fusion is the main production process of KK 
gravitons, that is, the A = ±2 states are mainly produced. Therefore, the azimuthal 
distribution is expected to be suppressed around $'^2 = 1^/2. However, for the higher 
graviton mass region, the contribution from qq annihilation, which enhance the events 
around $'^2 = '^/'^i becomes larger and gives rise to cancel the correlation each other. At the 
graviton mass M = 3.4 TeV the production rate of the two contributions is comparable [19], 
and hence the azimuthal distributions are enhanced around <I>'^2 = 7r/2 since the correlation 
in qq annihilation is stronger than that in gluon fusion; see ( ^.39 ) and ( 5.4C| ). 

In the future e^e~ (or photon) linear collider, on the other hand, only the A = ±1 (or 
A = ±2) state can be produced, and more precise studied to determine the spin-2 nature 
may be done in the clean environment. 
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6. Summary 



We have studied angular correlations of the two accompanying jets in Higgs boson and 
massive-graviton productions at hadron colliders, which include WBF and GF processes. 
We have also considered their decays into a vector-boson pair which subsequently decay 
into qq or gg. 

The amplitudes for the VBF subprocesses are given by the product of the two incoming 
current amplitudes and the off-shell VV X amplitude summed over the polarization of 
the t-channel intermediate vector-bosons (see eq. (2.7)), while the amplitudes for the X 
decays into four final states via a vector-boson pair are expressed as the product of the 
X VV amplitude and the two outgoing current amplitudes summed over the polarization 
of the s-channel intermediate vector-bosons (see eq. ( 2.15| )). Using the kinematical variables 
in fig. ^, we presented all the helicity amplitudes explicitly; tables |l] and § for the quark and 
gluon currents, and tables § and|5| for the VV — H/A and VV — > G processes, respectively. 
We also showed that our off-shell vector-boson current amplitudes reduce to the standard 
quark and gluon splitting amplitudes with appropriate gluon-polarization phases in the 
collinear limit (eqs. ( |4.18| ) and ( [4.20|) ). 

To validate our analyses, we demonstrated that the VBF amplitudes dominate the 
exact matrix elements not only for the WBF processes but also for all the GF processes 
when typical selection cuts to enhance the VBF events are applied, such as a large rapidity 
separation between two jets in (5.2). Furthermore, we found that the pxj slicing cut ( ^.^ ) 
is eff'ective to suppress the non-VBF diagrams especially for the graviton productions. 

By using the density matrix formalism in eq. (2.17) and our analytical amplitudes, we 
showed that nontrivial azimuthal angle correlations of the jets in the production and in the 
decay of massive spin-0 and -2 bosons are manifestly expressed in terms of the quantum 
interference among different helicity states of the intermediate vector-bosons; see (5.7) for 
the production and ( |5.28| ) for the decay. 

For the productions and the decays of Higgs bosons, our analytical arguments can de- 
scribe the previous studies on the angular correlations; for instance, the WBF gives the flat 
azimuthal distribution (|5.9|), while the GF produces the cos2((/)i — 02 ) distribution ( 5.1C| ). 
We also explicitly showed that the gluon currents, especially in the decay processes, re- 
duce the azimuthal correlations. For the massive-graviton case, we found the 0-dependent 
cos2((^i -|-</>2) correlations for both the production processes ( 5.18 ) and the decay processes 
( 5.34 ), which are proportional to sin^0. The correlations also depend on the final gravi- 
ton polarization in the productions and on the initial graviton polarization in the decays. 
Those correlations reflect the spin and the CP nature of the Higgs bosons and the massive 
gravitons, and may have a great potential to be observed at the LHC. 
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A. Wavefunction and vertices for a spin-2 particle 

The polarization tensor for a spin-2 particle in table ^, e^'^{p,X), is decomposed into po- 
larization vectors for a spin-1 particle as 

e^-(p,±2)=e^(p,±)6'^(p,±), 

6^-(p, ±1) = [e'^ip, ±) e-'ip, 0) + 6^(p, 0) e'^ip, ±)] , 
e'^-'ip, 0) = [e^^ip, +) e-'ip, -) + e^^ip, -) e'^ip, +) + 2 6^(p, 0) e'^ip, 0)] . (A.l) 

See more details in ref. [46]. 

The vertices for a massive graviton with two vector-bosons in table ^ are given by [56] 

f ^Ty (gi, 92) = {ml^ + qi • 92) C^'^'"^ + D^^^^P{q^,q2) + ^'£^^"'"^(91, 92), (A.2) 
where qi^2 and n-v are the momenta and mass of the vector bosons, and 

I?^"'°^(9l, 92) = g^'q^q^ - [g^^q'iq^ + g^'^qWi " ^"^f 92 + ^ ^)\ , (A.4) 
E^^'^f{q,,q2) = g^'iqU^^ + q^q^^ + qfq^,) - [g'^q^l + ^'^"'Zs^^^ + (/^ ^)\ ■ (A.5) 

The ^ term is the gauge-fixing term, which vanishes for massive vector-bosons in the unitary 
gauge 00). For massless vector-bosons we take ^ = 1 in the Feynman gauge. 

B. Relation between wavefunctions and d functions 

In this appendix, we demonstrate that two-to-two processes via s-channel spin-2 resonances 
can be factorized into the production part and its decay part, by using explicit spin-2 
wavefunctions. In the resonance (X) rest frame, p^ = (m, 0), we consider the processes 

ai{ki,ai) + a2{k2, (T2) X{p) a\{k[,a[) + a'^ik^, a'2), (B.l) 

where the oi momentum [ki) is taken along the positive z-axis and the a'^ momentum {k'^) 
is given by the scattering angle 0. 
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Before we consider the spin-2 resonance case, let us start with the well-known spin-1 
case. The numerator of the propagator for spin-1 particles, 



^ e^^ip,Xre''{p,X), 



(B.2) 



A=±1,0 
.2 _ ^2 



is a projector on the on-shell particle, p = m . Since the projector is an operator which 
satisfies = P, namely 



pM-(p) = P^^P{p)P^%p) 

= ^ e^ip, Xy eP{p, A) ep{p, A')* e^(p, A'). 



(B.3) 



Here, A is the helicity along the incoming ai momentum and A' is the helicity along the 
outgoing a'l momentum, that is, X = ai — a2 and A' = a[ — On the mass-shell, 
the propagator factor ( [B.SD can be factorized into the production part and the decay 
part; the wavefunction e^{p,X)* is used to calculate the production amplitudes, and the 
wavefunction e'^{p,X') is used to calculate the decay amplitudes. There is a connecting 
factor, eP{p, X) ep{p, X')* , which should be a scalar function depending on A, A', and the 
orientation angle between the two quantization axes. Using the explicit forms of the 
spin-1 polarization vectors (with the HELAS convention [38]): 



e^{p,X 



V2 



(0, ^1, -i, 0), 



e^(p,A = 0) = (0, 0, 0, 1), 



(B.4) 



and 



e''{p,X' = ±) = ^(0, =Fcos6', -i, ±sin( 
v2 

€^'{p,X' = 0) = (0, sin 61, 0, cos6i), 
the connecting factor can be expressed as the J = 1 d function. 



(B.5) 



eP{p,X)e,{p,X' 



^(l + cos0) --^sin^ i(l-cosi 



^(l-cos6') 4^ sin 6* i(l + cos( 



sin( 



-4,A'(^)- 



(B.6) 



The d function dictates the overlap of the angular momentum states between the initial 
0102 and the final a'^a2 state. Finally, the propagator factor (|B.3| ) can be rewritten as 



(B.7) 



P^''{p) = - Y.e>^{p.Xr d\^{e)e^{p,}:). 

A,A' 



Note that the sign comes from the orthogonal relation, e^(p, A*^'-*) e^(p, A^'^)* = — 1. 
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d].2^i = -d^2-i =F^(1 ± cos6')sm 



"±2,0 




Table 10: Explicit forms oi J — 2 d functions. 

Now, let us move to the spin-2 case. Similarly, the propagator factor for spin-2 particles 
is a projector on the mass-shell (p^ = m^): 



A=±2,±l,0 



(B. 



where e^^'^{p, A) is a spin-2 wavefunction. On the mass-shell, the summation over the helicity 
can be duplicated. 



= ^ e'^-'ip, Xy eP^p, A) e,^{p, A')* e°^(p, A'). 



(B.9) 



A A' 
By using the explicit forms of the spin-2 polarization tensors in eq. ( [A.l] ) and the polariza- 



tion vectors in eqs. (B^) and (|B.5| ), one finds that the overlap factor is the J = 2 d function 
(see table IC) as 

e''-(p,A)6p.(p,A')* = 4A'W- (B.IO) 

Therefore, the propagator factor for spin-2 particles ( p. 9] ) can be factorized into the pro- 
duction and the decay part as 



= ^ xy dlyie) e^^ip, A'). 



(B.ll) 



A, A' 



C. Angular distributions for G VV {f f){f f) 



In the heavy graviton mass limit (/?' = y 1 — Amy /M^ — > l) , where the only three ampli- 
tudes among the nine G — > VV amplitudes in table ^ are relevant, the differential distribu- 
tion of the massive-graviton decays, G — > VV — > (/i/3)(/2/4); for the initial polarization 
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A = ±2 along the z-axis is given by the expression 
d cos G d cos Oi d cos O2 d^i2 

~(1 + + C^){(1 + C?)(l + 4) + 4K1K2C1C2} + 2sS2g2 

=F 8c(l + c2){kiCi(1 + 4) + '^2(1 + c?)c2} 

+ 4s^SlS2{(l + C^)(kiK2 + C1C2) =F 2c(k2Ci + K1C2)} COS <I>12 

+ cos 2«>12 (C.l) 

with $12 = (/)! + (/>2, the abbreviations c = cos0, s = sin 6, Cj = cos^j and Sj = sin^j, and 
in eq. ( |5.38| ); see fig. ^ for the definition of the kinematical variables. Here, primes (') 
which indicate the decay variables in the text are omitted for simplicity. For the A = ±1 
state, the distribution reads 

d cos d cos 6*1 d cos 62 (i$i2 

-4(1 - c^){(l + c\){l + cl) + 4K1K2C1C2} + Ss^c^s\sl 

^ 16c(l - C^){KlCl{l + cl) + K2(l + C?)C2} 

- 16s^CSiS2[c(kiK2 + C1C2) =F («^2Cl + K1C2)] COS ^>12 

-4sS?s^cos2$i2, (C.2) 
where the same normalization is used as in eq. ( p.l[) . 
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